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Abstract 


The  work  presented  here  concerns  microwave  interaction 
with  plasmas  in  gas  discharges  and  semiconductors.  The  main  aim 
of  the  investigations  has  been  to  get  a  basic  understanding  of  some 
central  problems  in  this  field  that  have  applications  to  plasma  ampli¬ 
fiers. 

The  different  sub- reports  can  be  summarized  as  follows. 

I.  The  excitation  of  waves  in  a  waveguide  partially  filled 
with  a  cool  isotropic  plasma  is  studied.  Special  attention  is  paid  to 
cases  where  backward  waves  exist  and  the  influence  of  a  drift  motion 
of  the  plasma  is  considered. 

II.  A  circular  waveguide  is  filled  with  two  concentric  die¬ 
lectrics  and  with  an  infinitely  thin  metal  tube  between  the  dielectrics 
from  z  =  0  to  z  =  4  oo.  The  effect  of  the  discontinuity  at  z  =  0  on  a 
dominant  circularly  symmetric  TM-mode  incident  from  z  *  -  oo  .is  in¬ 
vestigated.  The  results  have  applications  to  plasma  waveguides. 

III.  Wave  propagation  in  plasma  waveguides  is  studied  theore¬ 
tically  and  experimentally.  Dispersion  curves  showing  the  existence  of 
backward-waves  are  given  for  the  surface-wave  dipole  modes.  The 
quasi- static  solutions  have  been  compared  to  the  exact  ones.  Experi¬ 
ments  have  been  performed  and  show  good  agreement  with  the  theory. 

IV.  Electron  plasma  resonance  absorption  is  observed  in  an 
n-type  crystal  of  Ge  at  35.  000  Mc/sec  and  liquid  hydrogen  temperature. 
The  crystal,  in  the  form  of  a  rod,  is  placed  across  a  rectangular  wave¬ 
guide  parallel  to  the  broad  walls.  No  static  magnetic  field  is  applied. 
Microwave  power  absorption  of  a  magnitude  and  temperature  dependence 
closely  agreeing  with  that  predicted  by  an  approximate  theory  developed 
by  the  author  is  measured  using  a  suitably  doped  Ge  specimen.  Obser¬ 
vation  of  plasma  resonance  absorption  gives  a  direct  method  to  deter¬ 
mine  the  charge  carrier  density  in  semiconductors. 

V.  It  is  shown  that  the  ordinary  theory,  which  describes  a 
magnetized  plasma  as  an  anisotropic  dielectric,  can  be  used  to  describe 
the  magneto-plasma  resonances  observed  in  semiconductors  and  metals. 


VI.  The  variation  of  plasma  density  and  electron-temperature 
along  the  positive  column  of  a  low-pressure  Hg-discharge  is  studied  by 
means  of  a  microwave  cavity  and  Langmuir  probes. 

VII.  The  noise- radiation  from  a  low-pressure  Hg-discharge 
inside  a  waveguide  has  been  measured.  The  scattering  and  absorption 
of  a  coherent  signal  has  also  beon  measured  and  compared  with  the 
noise-spectra. 

VIII.  Using  Kirchhoff 's  law  and  a  simple  hydrodynamic  model 
it  is  shown  that  thermally  excited  longitudinal  waves  give  rise  to  noise 
emission  in  the  microwave  region  from  cylindrical  gaseous  discharges. 
The  emitted  radiation  has  a  characteristic  spectrum  with  several  sharp 
peaks  at  the  resonance  frequencies  of  the  longitudinal  standing  waves. 
The  predicted  spectrum  seems  to  agree  well  with  experiments. 
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I .  EXCITATION  OF  BACKWARD  WAVE  MODES  IN  PLASMA  WAVEGUIDES 


Introduction 

It  is  known  that  in  waveguides  partially  filled  with  a  homogeneous  medium 
that  may  either  be  isotropic  or  anisotropic  there  exist  modes  with  negative 
dispersion  (backward  waves).  This  medium  can  be  an  ordinary  dielectric 
(CLARRICOATS  1960),  a  plasma  (TRIVELPIECE  1959),  (OLINER  1962),  or  a  ferrite 
(TRIVELPIECE  1961).  In  order  to  get  more  insight  into  the  nature  of  these  modes 
we  will  study  the  excitation  of  a  circularly  cylindrical  waveguide  containing  a 
homogeneous,  cool,  isotropic  plasma  column  (Fig.  1).  The  frequency  of  the 
exciting  signal  will  be  assumed  to  be  so  high  that  the  motion  of  the  ions  can  be 


neglected. 

We  will  first  consider  under  what  conditions  negative  dispersion  may  exist 
in  the  system  shown  in  Fig.  1  and  discuss  the  general  behaviour  of  the  modes  in 
this  waveguide.  The  modes  that  exist  in  this  system  are  in  general  hybrid  modes 
with  both  Ez-  and  Bz-components .  Only  modes  with  no  azimuthal  field  variations 
are  pure  TE-  or  TM-modes. 

The  power  flow,  P,  in  the  z-direction  through  a  cross-section  of  this 
waveguide  can  for  one  mode  be  written  in  the  following  form: 
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is  the  plasma  frequency. 


For  a  lossless,  passive,  uniform  waveguide  one  can  deduce  the  following 
relation  between  the  power  flow,  the  group  velocity  (jg)i  and  the  stored  energy 
per  unit  length  (W), 


2. 


provided  p  is  real  (CHORNEY  1962).  Backward  waves  occur  when  the  phase  velocity 
and  the  group  velocity  are  opposite  (3  <  0 ). 

Equation  (1)  together  with  Eq.  (2)  shows  that  backward  waves  cannot  exist 
for  TE-modes.  Nor  can  they  exist  for  TM-roodes  if  £  0.  By  performing  the  inte¬ 

gration  in  Eq.  (1)  it  can  be  shown  that  backward  waves  cannot  exist  for  TM-modes 
even  if  £  <  0.  Thus,  it  is  only  for  hybrid  modes  that  backward  waves  may  occur 
in  the  system  shown  in  Fig.  1. 

j<T 

The  case  where  the  fields  vay  as  e  in  the  (^-direction  will  now  be 
discussed.  The  characteristic  equation  for  this  case  has  been  studied  (ALHIEDS- 
SON  1963)  and  the  dispersion  curves  have  the  general  form  shown  in  Fig.  (2. 

The  dispersion  curves  I,  II,  ...  correspond  to  modes  which  are  very  similar 
to  those  in  the  empty  waveguide.  These  modes  are  infinite  in  number  and  have  a 
phase  velocity  larger  than  the  velocity  of  light.  The  waves  associated  with 
dispersion  curve  1-2-3  can  propagate  with  a  phase  velocity  smaller  than  the 
velocity  of  light  and  are  very  similar  to  waves  on  a  free  plasma  rod.  At  the 
points  and  Q^,  where  P'd/3 3  =  0,  the  dispersion  curve  1-2-3  will  be  connected 
with  the  dispersion  curves  1’,  2',  3'  and  1",  2”,  3”,  which  are  only  schematical¬ 
ly  indicated  in  fhe  figure.  These  dispersion  curves  correspond  to  complex  values 
of  the  propagation  constants  and  for  one  pair  of  branches,  e.g.  1',  2',  the  J3- 
values  are  complex  conjugate  for  a  givpn  value  of  w.  If  losses  are  introduced, 
branch  1  will  no  longer  be  connected  with  branches  2  and  2'  but  only  with  branch 
1'.  The  situation  will  be  the  same  for  point  Qg.  Thus  the  dispersion  curve  1-2-3 
breaks  up  into  three  different  dispersion  curves  corresponding  to  three  different 
modes,  which  for  convenience  will  be  called  the  plasma  modes.  In  case  of  zero 
losses,  the  fields  corresponding  to  the  dashed  dispersion  curves  carry  no  energy 
(CHORNEY  1962)  and  decrease  exponentially  from  the  source  of  excitation.  The 
number  of  plasma  modes  may  be  one,  two  or  three  (as  in  the  case  above).  Fig.  3 
shows  the  dispersion  curves  for  the  plasma  modes  for  three  different  values  of 
the  plasma  frequency  at  given  r./r_.  In  the  following  we  will  primarily  be 

i  2 

interested  in  the  excitation  of  the  propagating  plasma  modes. 


Excitation  by  a  Gap 

We  will  now  determine  the  fields  in  the  system  shown  in  Fig.  1  when  it  is 
excited  through  the  gap  and  the  tangential  E-fields  at  r  *  r2  for  |  z|  <  1  are 
given.  In  the  Apoendix  it  is  shown  that  for  z>  0  there  will  only  exist  waves 
with  |j^  >  0  and  for  z  <  0  only  waves  with  •<  0.  Thus,  for  modes  with  negative 
dispersion  the  phase  velocity  will  be  directed  towards  the  excitation  gap. 
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The  differential  equation  for  E  (and  B  )  is 
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The  Fourier  transforms  yz(r,p)  and  ?2(r,p)  are  defined  as 
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and  ¥z  (and  £z)  will  thus  satisfy  the  equation 
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The  Fourier  transforms  of  the  transverse  conponents  are  given  hy 
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In  the  plasma  region,  r  <  a,  the  solution  of  (5)  ie  of  the  forn 
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where  h^  =  (p2-  ep2)^  and  h^  =  (p2  -  p2)‘ 
The  boundary  conditions  are 
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Using  the  boundary  conditions  together  with  (8)  and  (9)  we  get  the  following 
equation  for  the  coefficients  A,  B....: 
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where  the  following  abbreviations  have  been  used 


*11  =  I1(h1r1) 


*21  *  I1(h2r1) 


*22  =  W2> 


sind  similarly  for  K^r  1^  and  K.^ . 

Solving  A  in  Eq.  (14)  and  using  the  inverse  Fourier  transform  on  Eq.  (8)  we 
obtain  the  following  expression  for  the  E^-field  inside  the  plasma  column: 
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A(p)  =  0  is  the  characteristic  equation  for  the  plasma  waveguide  (AGDUR  1962). 

The  other  field  components  can  be  obtained  in  an  analogous  way.  For  z  <  O 
this  integral  is  evaluated  by  a  contour  integration  along  the  contour  shown  in 
Fig.  4.  The  integral  along  the  semicircle  goes  to  zero  when  R  -»  50  and  the  value  of 
the  contour  integral  is  equal  to  the  integral  (15).  The  integrand  in  (15)  has  no 
branch  points  (the  branch  point  singularities  from  individual  terms  cancel)  but 
an  infinite  number  of  poles  which  are  indicated  in  Fig.  4.  One  can  distinguish 
between  three  different  types  of  poles:  1)  poles  on  the  imaginary  axis  which 
correspond  to  pure  cut-off  modes;  2)  poles  in  the  complex  plane  occurring  in 
complex  conjugate  pairs  which  correspond  to  the  dashed  curves  in  Fig.  2  and  which 
in  analogy  with  the  cut-off  modes  do  not  contribute  to  the  energy  flow;  3)  poles 
on  the  real  axis  which  correspond  to  propagating  modes.  If  losses  are  introduced, 
the  poles  on  the  real  axis  will  move  either  upwards  or  downwards  in  the  3-plane 
depending  on  whether  they  correspond  to  waves  with  <  0  or  —  >  0,  as  shown  in 
the  Appendix. 

As  an  example  we  study  the  excitation  of  the  waveguide  at  a  frequency  u  at 
which  only  plasma  modes  are  propagating  and  have  a  dispersion  characteristic  as 
shown  in  Fig.  5.  From  the  Appendix  and  the  figure  it  is  seen  that 
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In  the 

limiting  case  of  no 

losses  we 

therefore  have  to  1 

shown 

in  Fig.  6. 

The  value  of  the  integral  in  Eq.  (15)  is  then  obtained  from 
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assuming  that  all  poles  are  simple.  Double  poles  occur  for  points  on  the  disper¬ 
sion  curve  where 


,  3 A  , 

*3p  =  0  * 


On  the  basis  of  the  theory  developed  we  will  now  study  the  excitation  of 

r2 

plasma  waves  for  a  plasma  waveguide  where  —  =  5  and  v  r_/c  -  1.7.  The  dispersion 

P  v 

curve  for  this  case  is  given  in  Fig.  3.  As  seen  from  the  figure  it  is  possible  to 
excite  three  different  plasma  modes  in  this  case.  The  cut-off  frequencies  for  the 
waveguide  modes  are  well  above  the  cut-off  frequencies  for  the  plasma  modes  for 
these  parameter  values,  and  for  to/w  <  the  olasma  modes  will  be  the  only  pro- 

P  v  2 


pagating  modes  excited  in  the  waveguide. 

The  excitation  will  be  studied  when  the  electric  field  in  the  gap  is:  E_,  = 

=  b(z  -  0);  E  =0,  which  corresponds  to  a  unit  voltage  across  an  infinitely 

7o  1 

short  gap  and  gives  R  =  -=  and  Q  =  0.  Far  away  from  the  gap  the  cut-off  nodes 

V2 


have  a  negligible  amplitude  and  only  the  propagatin:  plasma  nodes  will  contribute 

to  the  fields.  Fig.  7  shows  the  E^-field  at  the  boundary  of  the  plasma  column  for 

the  plasma  modes  excited  by  the  unit  gap  voltage.  Ihe  power  flow  in  a  plasma 

mode  when  this  mode  has  a  given  value  of  the  E  -field  at  the  plasma  boundary  is 

z 

shown  for  the  three  plasma  modes  in  Fig.  8.  The  power  flow  has  been  computed  from 

Eq.  (1)  and  E  at  the  boundary  has  the  value  I„(h..r.).  By  combining  the  results 
z  ill 

shown  in  Fic.  7  and  Fig.  8,  Fig.  9  has  been  obtained  which  shows  the  power  flow  P 

in  the  different  modes  associated  by  the  above  excitation. 

The  asymtotes  in  these  figures  correspond  to  the  two  points  on  the  dispersion 

curves  that  seoarate  the  different  plasma  modes.  At  these  points  E  and  P  -*  °°. 

du  z 

P^  -»  0,  which  follows  from  Eq.  (2)  because  is  finite  while  mpr  -*  0.  As  pointed 
out  above  the  character  of  the  dispersion  curve  will  change  drastically  around 
these  points  when  losses  are  introduced  and  the  sane  will  hold  for  the  curves 
shown  in  the  last  figures.  When  p  -»  0,  -*  0  because  the  plasma  mode  becomes  a 

pure  TE-mode  in  the  limit  p  =  0.  Per  the  same  reason  P  -♦  *  when  3  -»  0.  P^  =0 
when  p  -*  oe  simply  because  the  given  E^-field  at  the  plasma  boundary  approaches 
infinity. 
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The  main  conclusion  that  can  be  drawn  from  the  results  in  Figs.  7  and  9  is 
that  the  power  and  electric  field  contributions  from  mode  C  are  in  general  very 
small  in  comparison  with  the  contributions  from  the  two  other  plasma  modes. 

Excitation  by  a  Given  Transverse  Field 


For  comparison  with  the  results  obtained  above  we  will  here  study  the  power 


flow  in  the  plasma  modes  for  the  sane  waveguide  as  above  (— ^  =  5;  -  1.7) 

r  1 

when  the  transverse  field  in  the  plane  z  =  O  is  given  by 
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where  is  the  first  root  of  J’(x)  =  0.  This  is  the  sane  field  as  the  transverse 

electric  field  for  the  TE  -mode  in  an  empty  waveguide  with  radius  r_ .  Expanding 

11  £ 

the  field  in  the  plasma  waveguide  for  z  =  0  in  the  different  nodes  and  allowing 

fto) 

only  modes  with  >  0  we  put 
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is  the  transverse  field  in  the  n:th  node.  Using  the  ortho- 
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we  obtain  from  Eq.(18): 
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The  power  flow  in  the  n:th  mode  will  then  be: 

_  x  5*  )psj2 


(21) 


Tn  z 


has  been  computed  for  the  plasma  nodes  treated  cn  page  and  the  results  are 
given  in  Fig.  10.  Also  in  this  case  the  power  transported  in  the  slowest  plasma 
mode  (C)  is  in  general  negligible  in  comparison  with  that  transported  in  the 
other  two  modes. 
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Influence  of  a  Drift  Velocity  on  Dispersion  Curves 

If  the  plasma  electrons  in  the  plasma  waveguide  treated  above  are  assumed  to 
have  a  drift  velocity  the  dispersion  curves  for  the  new  system  («',  p')  can  be 
obtained  by  the  following  transformation  of  the  ui, p-curves  given  above. 


vo 

(— ) 


(>j)  +  pvQ) 


(22) 


P* 


(23) 


is  the  drift  velocity  of  the  plasma  electrons  in  the  positive  z-direction.  The 
dispersion  curves  in  the  w, p- system  are  not  affected  by  a  motion  of  the  waveguide 
walls. 

Fig.  11  illustrates  the  influence  of  a  drift  velocity  of  the  electrons.  The 
dashed  curve  corresponds  to  =  0.  In  the  case  shown  in  this  figure  the  group 
velocity  in  the  u* , p' -system  is  always  positive  and  we  can  therefore  only  excite 
propagating  plasma  modes  in  the  positive  z-direction.  The  group  velocity  for  the 
plasma  modes  treated  above  may  be  very  small  and  even  a  very  snail  drift  velocity 
may  have  a  large  effect  on  the  dispersion  characteristic  of  the  system.  This  is 
shown  by  Fig.  12. 


APPENDIX 


The  waveguide  described  above  belongs  to  a  class  of  waveguides  which  have  a 
characteristic  equation 


F(u,p,jv)  =  0 


(24) 


where  F  is  an  even  function  of  6  and  real  for  real  values  of  u,  (5  and  jv.  v  is 
a  loss  parameter. 

If  “jt  P^  (both  real)  is  a  solution  to  Eq.  (24)  when  v  =  0  we  can  obtain  the 
following  two  solutions  when  small  losses  are  introduced  (v  *  Av) 


(I)  (D  =  a)  ;  p  =  p.  - 


cf 

^(j?) ) 


1  ,3F 


-  •  JAv 


^(jvV. 


(II)  w  =  »0  - 


p  =  Pi 


(25) 


(26) 
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The  second  solution  corresponds  to  a  cavity  making  damped  oscillations  and  we 
must  therefore  require 


Considering  the  behaviour  of  a  wave  e-J^z, excited  at  z  *  O,  when  z  it  is 

necessary  that  Irn(p)  <  O.  For  such  a  wave  we  therefore  have  the  following 
condition 


We  therefore  finally  obtain 


do) 

0P 


3F 

00) 


(29) 


This  means  that  only  those  waves  which  have  the  group  velocity  directed  outwards 
from  the  place  of  excitation  will  be  excited. 

This  result  has  been  proved  before  using  more  sophisticated  methods  (BRESLER, 
1959). 
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Metal  tube 


Fig.  1  Circularly  cylindrical  waveguide  with  central  plasma  column 
and  a  gap  in  the  waveguide  wall. 
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II.  A  DISCONTINUITY  PROBLEM  IN  A  CIRCULAR  WAVEGUIDE 
CONTAINING  TWO  DIELECTRICS 

Introduction 

In  the  study  and  application  of  plasma  waveguides  it  is 
important  to  know  how  energy  is  transferred  between  conventional 
waveguides  and  plasma  waveguides.  An  attempt  to  obtain  approximate 
solutions  of  this  problem  by  a  variational  method  has  been  made 
(KALMYKOVA  1961). 

In  this  paper  we  shall  investigate  a  waveguide  structure 
for  which  exact  solutions  can  be  obtained.  The  problem  is  formulated 
and  solved  for  a  waveguide  that,  for  values  of  the  dielectric  constant 
less  than  one,  represents  a  waveguide  containing  a  cool,  homogeneous 
plasma. 

Chapter  1,  Formulation  of  the  problem 

The  waveguide  system  which  we  shall  investigate  in  this 
paper  is  shown  in  Fig.  1.  The  part  from  i  =  Oto  z  =  -  ®  consists  of 
a  circularly  cylindrical  metal  tube  containing  two  concentric,  lossless 
dielectrics.  The  dielectric  constant  is  in  the  region  r  <  a  and 
in  a  <  r  <  b.  In  the  following  this  waveguide  will  be  called  waveguide  1 . 
The  waveguide  system  from  z  =  0  to  z  =  +  ®  is  identical  with  the  first 
part  except  for  an  infinitely  thin  metal  tube  between  the  two  dielectrics. 
This  part  thus  consists  of  two  waveguides:  a  coaxial  line  (waveguide  2) 
filled  with  dielectric  (  e^),  and  a  circularly  cylindrical  waveguide 
(waveguide  3)  with  radius  a  and  filled  with  dielectric  (e  ^). 

We  shall  consider  the  effect  of  the  discontinuity  at  z  »  0 
when  the  dominant,  circularly  symmetric  TM-mode  in  waveguide  1  is 
incident  from  z  =  .  Because  of  the  discontinuity  the  incident  wave 

will  generally  excite  all  circularly  symmetric  TM- modes  in  the  three 
waveguides,  including  the  TEM-mode  in  waveguide  2.  For  each  mode. 
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only  the  wave  with  the  group  velocity  directed  away  from  the  discon¬ 
tinuity  will  be  excited. 

By  means  of  a  method  based  on  the  WIENER-HOPF  technique 
* 

(see  e.  g.  JONES  1952,  HALLEN  1953)  we  obtain  the  solution  to 
MAXWELL*  s  equations  that  gives  the  scattered  field  due  to  the  dis¬ 
continuity,  The  general  solution  has  such  a  complex  mathematical 
form  that  it  is  of  little  practical  value.  Simpler  solutions  can,  however, 
be  found  for  certain  conditions.  Thus  we  shall  restrict  the  values  of 
the  parameters  e^,  e^,  a,  b,  and  the  frequency  of  the  incident  wave, 

(u,  so  that  a)  in  waveguide  1  only  one  mode  can  propagate  while  the 
other  modes  are  cut  off,  b)  in  waveguide  2  all  modes  except  the  TEM- 
mode  are  cut  off,  c)  in  waveguide  3  all  modes  are  cut  off.  The  fields 
of  the  cut-off  modes  decrease  exponentially  away  from  the  disconti¬ 
nuity.  For  large  z-values  they  become  negligible  compared  with  the 
fields  of  the  propagating  modes.  As  the  system  is  lossless,  we  now 
need  only  three  parameters  to  completely  characterize  the  fields 
far  from  the  discontinuity,  e.g.  the  amplitude  and  phase  of  the  reflec¬ 
ted  dominant  mode  in  waveguide  1  and  the  phase  of  the  transmitted 
TEM-mode  in  waveguide  2. 
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Chapter  2.  Some  characteristics  of  the  relevant  mode* 

In  this  chapter  we  mainly  consider  how  the  conditions  a,  b, 
and  c  of  Ch.  1  can  be  fulfilled. 

For  TMjj^-modes  in  waveguide  3  with  a  time-  and  z-depen- 
dence  e- P3mz)  we  have  the  following  dispersion  relation: 


or 


)<*/‘iko  -  “L>  = 0 


Pm\2 


(0 

(2) 


where  ^3m  wave  number  for  the  TMq^ -mode 

Jq  is  the  zero  order  BESSEL  function 

k_=“ 

0  c 

c  is  the  velocity  of  light 

p  is  the  m  root  of  Jn(p)  =  0 
*m  0'r' 

All  modes  are  cut  off  in  this  waveguide  if 
pl  C 

03  <  -  or  e  .  <0 


Vs  t  a 


1 


In  waveguide  2  the  TEM-mode  has  the  wave  number 

03  Vs  ,  - 


(3) 


(4) 


and  propagates  if  >  0.  If  the  wave  number  for  the  TMgm-mode  in 
this  waveguide  is  denoted  by  P^  m-K  We  ^ave: 

(5) 


«2  ,2  2 
^2,  m+1  e2  0  ” 


th 


where  g  is  the  m  root  of  the  equation: 
m 


J0(as)  Y0(bg)  -  YQ(ag)  J 0(bg)  =  0 


where  Yq  is  the  zero  order  NEUMANN  function. 


(6) 


26. 


All  TMgm- modes  are  cut-off  if 


The  dispersion  equation  for  TMq^- modes  in  waveguide  1 
is  (see  e.  g.  AGDUR  &  ENANDER  1961): 

£2J0(ahi)  [j|j(ah2)Y0(bh2)  -  YJj(ah2)J0(bh2)]  - 

c  h  <8> 

-  iq-1  Ji<abl>  lJ0<ah2>Y0<bh2'  -  Yolah2>I0<bh2>l  ■" 

where  hj  =  (c  ^  -  Pj)  2 

The  wave  numbers,  Plm»  in  this  waveguide  are  those  values  of  Pj  that 
satisfy  Eq.  (8).  In  order  to  fulfil  condition  a  we  must  assume  that 
only  one  of  these,  say  P^,  is  real. 

The  wave  numbers  for  the  cut-off  modes  in  waveguide  1  may 
be  complex^  as  distinguished  from  the  wave  numbers  for  the  cut-off 
modes  in  waveguides  2  and  3  which  are  purely  imaginary.  If  P^n  is  a 
complex  wave  number  that  satisfies  Eq.  (8)  then  the  complex  conjugate 
value,  P*n>  also  satisfies  this  equation.  This  can  be  proved  by  means 
of  the  principle  of  reflection  for  analytic  functions. 


1) 


The  proof  by  ADLER  (ADLER  1952),  that  the  wave 
numbers  of  a  waveguide  of  this  kind  cannot  be  complex, 
is  not  correct. 
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Figure  2  shows  dispersion  curves  for  TMQm-modes  in 
waveguide  1  for  a  case  where  the  inner  dielectric  is  a  plasma  and  the 
outer  is  vacuum.  In  this  diagram  complex  branches  of  the  dispersion 
curves  (not  shown  in  the  figure)  start  out  from  points  where  the  curves 
have  a  vertical  slope,  e.  g.  and  A^.  The  complex  branches  are  in 
this  case  connected  to  branches  where  the  wave  number  is  purely 
imaginary  as  distinguished  from  the  complex  branches  found  for  dipole 
modes  in  similar  waveguides  (AGDUR  1963,  OLINER  1962,  AGDUR  & 
ENANDER  1963). 

In  the  case  shown  in  Fig.  2  the  dominant  mode  has  a  zero- 
cut-off  frequency  and  a  phase  velocity  smaller  them  the  velocity  of  light. 
This  mode  resembles  the  surface-wave  mode  on  a  plasma  column  and 
will  always  exist  when  Sj  <  -c  ^  and  e  2  >  0.  When  both  e  ^  and  t  2 
are  positive,  the  dominant  mode  has  a  cut-off  frequency  different  from 
zero  and  the  dispersion  diagram  resembles  that  of  an  ordinary  circular 
waveguide  (BANOS  1951). 
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Chapter  3.  Analysis 

Assuming  a  time  dependence  e"*01,  we  find  the  following 
differential  equation  for  the  z- component  of  the  total  electric  field  in 
the  coordinate  system  shown  in  Fig.  1. 


82E 


8r 


zt 

7T 


,  8E  „  82E 

+  -  — +  ek2E  + - 

r  ST  0  zt 


=  0 


(9) 


,  ,2  ix> 

where  kg  =  — ^ 

c, 

and  t  =  £j  when  r  <  a 

e  =  c  ^  when  a  <  r  <  b 


The  boundary  conditions  are: 


Ezt  =  0  when  r  =  b,  -  oo  <  z  <  +oe  (10) 

Ezt  =  0  when  r  =  a,  0<z<  +  oo  ( 1 1 ) 

Ezt  continuous  when  r  -  a,  -ao  •<  z  <  0  (12) 

The  edge  conditions  at  the  discontinuity  give  (BOUWKAMP  1950) 

E  .  •*  C.  z  ^when  r  =  a  and  z  -*  -0  (13) 

zt  1  '  ' 

We  now  divide  the  total  electric  and  magnetic  fields  into 

two  parts,  one  representing  the  incident  field  and  one  the  scattered 

field,  so  that 


E  *  E  .  +  E 
Zt  Zl  z 


(14) 


In  the  whole  region  -oo  <  z  <  +oo,  r  <  a,  the  field  of  the  incident 


TMq  j-mode  is  taken  as 


Ezi-J0<rhll>e' 


jpil* 


(15) 


where  is  the  propagation  constant  for  the  dominant  mode  in  wave¬ 
guide  1  and 

hll  =  <‘lk0  "  pii^2 
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The  scattered  field,  E  ,  shall  satisfy  the  differential  equation  (9)  and 

z 

the  conditions  (10),  (12)  and  (13)  but  not  (11).  Instead  of  (11),  we  require 

jP..* 

=  -E^  =  -  J^rhjjJe  when  r  =  a,  0  <  z  <  +oo  (16) 

We  shall  now  introduce  small  losses  in  the  dielectric  by 
letting  c  j  and  e  ^  be  complex.  The  two  propagating  modes  are  then 
slightly  attenuated.  The  losses,  which  may  be  arbitrarily  small,  are 
introduced  to  facilitate  the  mathematical  procedure;  we  shall  eventually 
let  them  tend  to  zero. 

The  wave  numbers  for  the  two  propagating  modes,  P}1  and 
P21,  will  now  have  small  imaginary  parts  (smaller  than  the  imaginary 
part  of  any  wave  number  for  the  cut-off  modes): 

Pll  =P*ll+jPi’i 

11  11  11  (1?) 

P21  =  P2l  +  j(321 

where  Pj  j ,  P^j,  P j’j ,  and  p^  are  real  positive  numbers  (see  Fig.  3). 
Choose  6  equal  to  the  smallest  one  of  Pj’j  and  P^’j.  Then 

)E  I  <  C  .•  when  z  -*•  -oo  (18) 

Z  1 


iEzl  <  C2  *  e" 


when  z  -*•  +oo,  r  >  a 


<  •  e~^  ^  Z  I  when  z  -»  +oo,  r  <  a 


Introduce  the  FOURIER  transform  of  E 


vb  =  — L.t..  f  £  e^Q 

(2t)  1  Z 


jaz  . 
eJ  dz 


Transformation  of  Eqs.  (9),  (10),  (12),  and  (15)  yields 

+  7  S  +  (*ko  - a2) +  =  0 

ar 


=  0  at  r  =  b 
4;  continuous  at  r  =  a 
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The  solution  of  Eq.  (22)  that  satisfies  the  boundary  condi¬ 
tion  (1 5)  can  be  written 


4»  =  A(a)  In(7.r) 


for  0  <  r  <  a 


4-  =  B(a)  [I0(?2r)  KQ(72b)  -  K^r)  I0(72b)] 


for  a  <  r  <  b 


(25a) 


(25b) 


where  Ig  and  Kg  are  the  modified  zero-order  BESSEL  functions  and 

,  2  .  2v| 

71  =  (a  -CjV 

2  2  i  (26) 

72  =  («x  -  e2 k2)2 

7^  and  72  are  made  single-valued  through  suitable  branch  cuts  in  the 
complex  a -plane. 

For  circularly  symmetric  TM-modes  we  have 
3E  8ZH 

~  “  j®*1 0Hcp  +  (27) 

H  shall  be  continuous  at  r  =  a  for  -oo  <  z  <  0 
T  - 

and  satisfy  the  edge  condition 

i 

H„  -*•  C,z2  when  r  =  a  and  z  -*  0  (28) 

<p  2 

Denote  the  FOURIER  transform  of  H  by  4.  Transformation  of  Eq.  (27) 


gives: 


’  0  ddi  . 


where  7  =  7j  for  r  <  a  and  7  =  72  for  r  >  a 
Using  Eqs.  (25a)  and  (25b)  we  get: 


j<u*  .'o 

$  = - - -  A(a)  I^(7jf)  for  0  <  s  <  a 

j“«  2C  0 

®  = - —  B(a)[IJ,(72r)  KQ(72b)  - 

"  K0^72r^  V72b)J  for  a  <  r  <  b 
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Divide  tjj  and  9  into  two  parts 


where 


4*  =  4<+  + 

CO 

+t « -i-r  ? 
*  (2')'  i 


9  r  4>+  +  $_ 


E*  ejaZ,  i|i  = — i-t -  \E  ejazdz  (33) 
z  dzj  -  (2*)*  z 


and  4>+  and  #  are  similarly  defined. 

From  the  behaviour  of  the  fields  at  z  =  Joo,  Eqs.  (1^,  (19), 
and  (20),  we  conclude  that  and  4>+  are  analytic  functions  in  the  half¬ 
plane  Im(a )  >  -6,  and  that  ij;  and  are  analytic  for-Im(a)  <  6. 

At  r  =  a,  ijj  ,  and  9  shall  be  continuous,  whereas  9+ 
is  discontinuous.  Let  the  values  at  r  =  a  be  i|>+(a),  4'_(a)»  *_(a)»  ®+(a+0)» 
and  4+(a-  0).  ij/+(a)  is  obtained  from  Eqs.  (16)  and  33): 

i  Vhl  1  ^ 

(a.)  =  -  — i—tr  .  g — j-  (34) 

+  (2t)2  va  +  ^l  i ' 

From  Eqs.  (17),  (18),  (22),  (23),  (24),  and  (25),  we  now  derive 

i  J0(hlla) 

(35) 

=  B  (I0(72a)  KQ(72b)  -  K0(72a)  I0(7zb)] 

*+(a+0)  +  4>_(a)  =  -  B  -J  °  [l^a)  K^b)  - 

2  (36) 

-  Kb^2a)  I0^2b)] 

♦+(a- 0)  +  #_(a)  = - ~~  A  I^7ja)  (37) 

By  subtracting  Eq.  (37)  from  Eq.  (36)  and  eliminating  A  and  B  from 
the  resulting  equation  by  means  of  Eq.  (35)  we  get 


A+(a)  L(a)  - 


j  Jo<hna> 


where  we  have  put: 
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L  =  { 1/j <a«  Q]  •  (iQ^a)  Il0(72a)  K0^72b^  " 

"  K0^2a^  ^I0^2a^  K0^2b^  " 

"  K0^2a^  M72b^  1  ”  (€  \  ? 2^0  *0^1 a^  ^0^2a^  ^0^2b^ 

-  K0^2a> 

and  4>+(a-  0)  -  ®+(a  +  0)  =  A+(a )  which,  incidentally,  is  the  FOURIER 
transform  of  the  surface-current  density  in  the  inner  conductor. 

In  Eq.  (38),  which  is  valid  only  in  the  strip  -6  <  Im(o)  <  6 
(see  Fig.  3),  both  A^a)  and  (a)  are  unknown  functions  of  a ,  Never¬ 
theless,  they  can  be  evaluated  from  this  equation  by  means  of 
LIOUVILLE'S  theorem.  To  do  this  we  shall  rewrite  Eq.  (38)  in  such 
a  form  that  one  side  of  the  equation  is  regular  and  bounded  in  the  half¬ 
plane  Im(a )  >  -6  and  the  other  side  is  regular  and  bounded  in  the 
half- plane  Im(a)  <  6. 


We  introduce 

(a-P.i)  (a+Pn) 

«(->  “CH^iTT  <« 

and  assume  that  M(a)  can  be  factored  in  the  following  way: 

M(a)  =  M+(a)  M_(a)  (41 

wher«3  M+  is  regular  and  non- zero  for  lm(a)  >  -6  and  M  is  regular 

i 

and  non-zero  for  Im(a)  <  6  and 

M+(a)  ~  (a  2)  when  |  a  j  -*•<»,  Im(a)  >  -6 


M  (a )  ~  (a  ”  2  )  when  j  a  |  -*  oo,  Im(o )  <  6 
Using  Eq.  (40)  we  write  Eq.  (38)  in  the  following  form 


A  (a)  L  (o)  + — i-t 
*  +  (2*) 1 


Vhna> 


-'^1 1M“  1J 
j ,  yhna>  y  i 


TTWJ 


7T  ITF7 “  TT-PjTT 
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a+p«l  a-0n 

where  L+(a)  =  — ^ —  M+(a)  and  L_(a)=^--p —  M_(a) 
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The  left-hand  side  of  this  equation  is  a  regular  function  in  the  half-plane 
Im(a)  >  6  and  the  right-hand  side  is  regular  in  the  half-plane  lm(a)  <  6. 
The  two  half-planes  have  a  common  strip  and  the  two  sides  of  Eq.  (34) 
together  define  a  function  that  is  regular  in  the  whole  a-plane.  Further, 
we  find  from  the  edge  conditions  (13)  and  (28)  and  the  definitions  of  tjj  , 
Eq.  (33),  and  A+  that 

A+  ~  when  |a  j  -*•<»,  Im(a)  >  0 

-1/2  ,  ,  f*4) 

~  (a )  '  when  |  a  j  —■  oo ,  Im(a )  <  0 

Combining  this  result  with  the  assumed  behaviour  of  M+  and  M  at 
infinity  we  find  that  both  sides  of  Eq.  (43)  tend  to  zero  in  their  respec¬ 
tive  half-planes  when  |a|  oo.  According  to  LIOUVIL.LE  'S  theorem 
the  integral  function  defined  by  the  two  sides  of  Eq.  (43)  must  then  be 
identically  zero.  Thus,  putting  the  right-hand  side  equal  to  zero  we 


(1,  i  Jo<hna>  L»  , 

_(a)  -  —  i  j—rn — r  tth — r  J 


and  from  Eq.  (35) 

AI0(7la)  =  B[l0(r2a)  KQ(72b)  - 

(46) 

,  i  Vbiia^  L 

-  K0(72*)  *o<7Zb>  j  =  '  —f  (a+Fn)  L_(-Pn) 

The  waveguide  fields  can  then  be  obtained  by  using  the  FOURIER  inver¬ 
sion  theorem  on  Eqs.  (25). 


V 
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Before  doing  this,  we  shall  construct  functions  M+  and  M 
having  the  require^  properties.  Since 

Kq(z)  =  -  IQ(z)  ln(z)  +  f(z)  (47) 

where  f(z)  is  an  even,  integral  function,  the  ln(z)  terms  in  L,  (Eq.  (39)), 

cancel;  thus  L  has  no  branch  points.  The  numerator  in  the  definition 

of  L  is  an  integral  function  of  a  which  is  zero  if  and  only  if 

+ 


a  =  “  P-3 


/a  or  ci  “  •  k> 

2m  3m 

where  an<*  ^3m  are  t*le  wave  numbers  for  waveguides  E  and  3.  The 

denominator  is  also  an  integral  function  of  a  and  is  zero  if  and  only  if 

a  =  -  Pjm»  where  are  the  wave  numbers  for  waveguide  1.  M(a) 

has  the  same  properties  as  L(a)  except  that  M(a)  does  neither  have  zeros 

at  a  =  -  nor  F'les  at  a  =  -  We  now  define  M+(a)  and  M  (a)  by 

oo  +  jd 


Po 


and 


dz 

z  -  a 

(48) 

-ao  +  jd 

to  +  je 

dz 

z  -  a 

(49) 

-oo  +  je 


where  -6  <  d  <  e  <  +6 

The  use  of  CAUCHY  principle  values,  indicated  by  P  in  front  of  the 
integral  sign,  is  necessary  because  M(u)~  |a  |  when  a  -*  oo.  in  the 
strip,  so  that  the  ordinary  integrals  are  divergent* 

Defined  in  this  way,  M+  will  be  regular  and  non- zero  in  the 
half-plane  Im(a )  >  d.  Similarly  M  will  be  regular  and  non-zero  for 
Im(a)  <  e.  Further,  from  CAUCHY'S  integral  formula  we  find  M+M  = 
M,  when  -6  <  Im(a)  <  6. 


35. 


To  investigate  the  behaviour  of  M+(a  )  as  J  a  ]  -»  oo  in  the 
upper  half- plane,  we  write 
oo+jd 

P  \  ^  dz  =  I.  +  I_  +  I, 


where 


J  z  -  a 
-oo+jd 

D+id 


12  3 


u-Md 

1  =  J  TT?  dz  =  °^)  when  1°  I  - 


00 


-D+jd 


D 

»  ln  M(x  +  id)  i  M(-x  +  jd) 

I  =  \  x  + _ x 

3  J  x  +  jd  -  a  -x  +  jd  -  a 


It  can  be  shown  that 

12  =- jin  (a  |  +0(1) 

13  =  0(1) 

when  |  a |  -*•  oo;  Im  a  >  k,  where  d  <  k  <  0 

(50) 

when  |  a  |  -*  oo;  Im  a  <  1  where  e  > 

>  1  >  0 

so  that  M+  and  M  ,  defined  by  Eqs.  (48)  and  (49),  have  the  properties 
assumed  earlier. 

The  fields  in  the  waveguide  system  can  be  obtained  from 

Eqs.  (25)  and  (46)  by  means  of  the  FOURIER  inversion  formula.  As 

pointed  out  in  the  introduction,  it  is  usually  only  the  propagating  modes 

that  are  of  interest.  The  scattered  E  -field  at  any  point  in  the  region 

z 

r  <  a  is  given  by 


Therefore 


M+  (a)  =  0(  [a  |2  ) 


Similarly  we  find 

1  1 


WIT)  =°<laP) 
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oo  +  jf 

E  =— L*  f  A(a)I  (7  r)e"jaZda  (51) 

(2»)2 

'  '  -00  +  jf 

where  -6  <  f  <  6 

For  z  <  0  we  evaluate  the  integral  in  Eq.  (51)  by  closing  the  integration 
path  in  the  upper  half-plane  and  using  the  calculus  of  residues. 

From  Eqs.  (43),  (36),  and  (32)  we  then  obtain: 

i  C  in  °  *  Pi 


^  j  Jo(hila)  zpll  f  L(a)  a_P21 

Ez  =  -  *  rpjj)  *  F-'+^zi  j 

V7lr)  e-jQZ 

*  VvJ  "57^u 


(52) 


da 


The  integrand  in  (52)  has  poles  at  and  only  at  a  =  and  the  residues 

at  these  poles  give  the  amplitudes  of  the  reflected  modes  in  waveguide  1. 
The  reflected  wave  in  the  dominant  mode  is  thus: 


^  _PU_P21  M(Pn)  „  ^ 

Ezd"  Fu  +021  M_i-Fn}M+(Fn;  Vhnr>e 

For  the  reflection  coefficient  R,  defined  as 

R  =  (^W 

Z1 


-jpllz 


(53) 


we  get  from  Eq.  (53) 


R  = 


P11  ‘  P21 


21 


M  (Pu) 

3TO 


(54) 


From  Eqs.  (48)  and  (49)  we  see  that  when  the  losses  tend  to  zero 


M_(PU)  =  M+(Pn) 


so  that 


|R !  =1 


P11  ”  P21 

*77^ 


x  11 ' 

^11  +  *21 


(55) 


where  X  ^  ^  and  X£^  are  the  wavelengths  for  the  two  propagating  modes. 
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For  z  >  0  we  can  obtain  the  fields  in  an  analogous  way  by 
closing  the  integration  contour  in  the  lower  hal£»plane.  The  Enfield 
will  contain  two  propagating  modes:  one  with  the  wave  number  Pjj, 
which  exactly  cancels  the  incident  wave;  the  other  with  the  wave~janmber 
is  the  transmitted  TEM-mode.  The  phase  of  the  TEM-mode 
together  with  the  complex  reflection  coefficient  (Eq.  (54))  completely 
characterize  the  fields  far  away  from  the  discontinuity. 

The  simple  formula  for  the  magnitude  of  the  reflection 
coefficient  Eq.  (55)  agrees  with  that  obtained  for  the.  special  case  when 
s  j  *  t  2  (MARCUWITZ  1951)  with  the  preliminary-  results  of  an  experi¬ 
mental  investigation  that  is  being  performed  at  our  laboratory. 


33. 

Bibliography 

1.  ADLER,  R.B.:  Waves  on  Inhomogeneous  Cylindrical 
Structures.  Proc,  I. R. E.  40,  (1952),  pp.  339-348. 

2.  AGDUR,  B.:  Notes  on  Waves  in  Plasmas.  Symposium  on 
Electromagnetic  Theory  and  Antennas,  Copenhagen,  June 
1962.  Pergamon  Press  1963,  pp.  301-303. 

3.  AGDUR,  B.  &  ENANDER,  B. :  Resonances  of  a  Microwave 
Cavity  Partially  Filled  with  a  Plasma.  J.  Appl.  Phys.  33_ 
(1962),  pp.  575-581. 

4.  AGDUR,  B.  &  ENANDER,  B. :  Excitation  of  Plasma  Wave¬ 
guides  with  Backward  Wave  Modes.  Ericsson  Tech.  19 
(1963):2,  pp.  251-268. 

5.  BANOS,  A.,  SAXON,  D.S.  &  GRUEN,  H. :  Propagation 
Characteristics  in  a  Coaxial  Structure  with  Two  Dielectrics. 

J.  Appl.  Phys.  22  (1951),  pp.  117-123. 

6.  BOUWKAMP,  C.  J, :  Diffraction  by  Small  Circular  Disks  and 
Holes.  Philips  Res.  Rep.  5  (1950),  p.  401, 

7.  HALLEN,  E. :  Elektricitetslara.  Stockholm  1953  . 

English  translation:  Electromagnetic  Theory.  London  1962  . 

8.  JONES,  D.  S. :  A  Simplifying  Technique  in  the  Solution  of  a 
Class  of  Diffraction  Problems.  Quart.  J.  Math.  (2)_3  (1952), 
pp.  189-196. 

9.  KALMYKOVA,  S.S. :  On  the  Problem  of  Matching  Plasma 
Waveguides.  Zh.  Tekh.  Fiz.  31_  (1961),  pp.  1374-1378. 
English  translation  in  Soviet  Physics- Tech.  Phys.  ^  (1962), 
pp.  1001-1004. 

10.  MARCUVITZ,  N. :  Waveguide  Handbook.  Radiation 
Laboratory  Series  10  .  New  York  1951, 

11.  OLINER,  A.  A.  &  TAMIR,  T. :  Backward  Waves  on  Isotropic 
Plasma  Slabs.  J.  Appl.  Phys.  33  (1962),  p.  231. 


\\\ 


V  ! 


. :  / . 


Fig.  2.  t»  -  P  diagram  for  a  circularly  cylindrical  waveguide 

containing  a  central  plasma  column,  b/a  ■=  4;  b/c  *  1 
(<o  =  plasma  frequency).  Dotted  curves  represent  real 
values  of  p.  Solid  curves  represent  imaginary  values  of  p. 
For  comparison  the  o-  P  diagram  for  an  empty  waveguide, 
radius  =  b,  is  shown  by  the  dctted  curves. 
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III.  INVESTIGATION  OP  BACDfARfV-WAVB  MOOES  IN  PLASMA  WAVEGUIDES 


Introduction 

In  recent  years,  swe  authors  (TRIVELPIECE,  1959  and  1961,  CLARRICOATS,  1960, 
NAPOLI  and  SWARTZ,  1963,  OLINER,  1959,  AGDUR  and  ENANDER,  1962)  have  studied 
wave  propagation  in  inhomogeneous  waveguides.  Using  ferrites,  plasmas  and  ordinary 
dielectrics,  it  is  possible  to  achieve  propagating  backward-wave  modes  in  non¬ 
periodic  waveguides.  This  paper  describes  the  surface-wave  dipole  modes  in  the 
plasma  waveguide  systems  shown  in  Pigs. 1  a-c. 

Some  experiments  have  been  performed  which  generally  show  good  agreement 
with  the  theory. 


Theory 

I.  General  Remarks. 

Wave  propagation  is  studied  in  three  different  waveguide  systems,  infinite 
in  the  direction  of  propagation  (2-direction),  with  the  cross-sections  shown  in 
Figs.  la-c. 

The  plasma,  which  is  taken  to  be  cool  and  lossless,  is  characterised  by  a 

2  2 

dielectric  constant  e  *  1  -  <0  /«  where  w  is  the  olasma  frequency.  Further, 

P  p  P 

the  plasma  is  assumed  to  be  homogeneous,  i.e.  no  radial  or  axial  density  varia¬ 
tions  are  taken  into  account.  The  signal  frequency  is  assumed  to  be  so  high  that 
the  motion  of  the  ions  can  be  neglected. 

The  electromagnetic  fields  are  assumed  to  be 


ai  Fi(r) 


\  Gi(r)  }• 


jn<r  j  (wt-pz) 


•  e 


(1) 


where  a.  and  b.  are  constants,  F.(r)  and  G. (r)  are  Bessel  functions,  8  is  the 
11  xi 

propagation  constant  and  n  characterizes  the  azimuthal  variation. 

The  dispersion  equations  are  derived  by  matching  the  electromagnetic  fields 
in  the  different  regions  using  the  ordinary  boundary  conditions.  This  leads  to  a 
linear,  homogeneous  system  of  equations  in  the  unde ter aimed  amplitude  coefficients 
Non-trivial  solutions  exist  when  the  system  determinant  equals  zero.  As  soon  as 
n  /  0  there  exist  no  pure  TE-  or  TM-oodes. 


II .  Plasma  column  surrounded  by  a  homogeneous  dielectric  enclosed  in  a  metal  tube. 

The  system  shown  in  Fig.  1  a  consists  of  a  central  plasma  column  surrounded 
by  a  homogeneous  dielectric  (e  )  and  a  metal  tube  (0  *  ■*).  The  dispersion 
equation  for  this  system  has  been  derived  by  several  authors.  It  is  repeated  here 


as  given  by  AGDOR  and  ENANDER,  (1962) 


[CP  J11k1  <J22Y21  “  V22J21)  ~  £r(J22Y21  '  J21Y22}]X 

j*  t 

*  [j11k1  (J22Y21  "  Y22J21)  "  (J22Y21  "  Y22J21>]  * 

2  k  2 

n  r'  [“£  "  ^i]  [r^]  (J22y2l“  Y22J21)(J22Y21  "  J21Y22) 


2 

where  J,  Y,  J  an(j  y’  represent  Bessel  functions  of  the  first  and  second  kind  and 
their  derivatives.  The  following  abbreviations  have  been  used: 

J11  =  -V^l*  J21  =  Jn(k2r1)  J22  ~  3rS*Zr2) 

and  similarly  for  Y,  J’  and  Y’ .  Further 

k.  *  k  (e  -  t2)^ 

1  op 

k  *  k  (C  ~ 

2  or 

k  =  <o/c 
o 

T  a  P/ko 

The  field  components,  omitting  the  factor  are 
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E2  =  CJn(k1r)  *  ***  r  -  r1 
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where  the  relations  be  tween  the  unknown  amplitude  coefficients  are  determined  by 
the  boundary  conditions. 

To  find  the  expressions  for  the  dipole  nodes  we  put  n  *  1  1  in  eqs.  (2)  and 

(3). 

Pig.  2  shows  the  principal  behaviour  of  the  dispersion  curves  obtained  fro* 
eq.  (2).  The  velocity  of  light  in  the  surrounding  medium  is  v  *  cAfc^.  This 
gives  the  boundary  between  fast  waves  (v^  >■  v)  and  slow  waves  (v^.  <  v)  which  is 
indicated  by  a  thin  solid  line.  The  figure  shows  two  different  types  of  curves 
E  -  P  ...  and  A  -  S  -  C  -  D,  indicating  that  there  exist  two  different  kinds  of 
modes  that  can  propagate  along  the  waveguide.  The  first-»entioned  curves  - 
infinite  in  number  -  are  fast  modes  related  to  those  existing  in  the  empty  wave¬ 
guide.  In  the  following  we  shall  mainly  discuss  the  mode  represented  by  A  -  B  - 
C  -  D.  This  is  a  surface-wave  mode  (plasm  mode),  it  is  closely  related  to  the 
modes  on  a'  free  plasma  column  (discussed  in  section  IV). 

Points  on  the  (u, p)-curve  where  5w/3p  »  0  are  branch  points  which  means 

that  the  dispersion  curve,  F(td,p)  =  O,  enters  the  complex  p-olane  along  two 

x) 

complex  conjugate  branches.  The  branches  (A*,  B‘  ...)  are  dram  schematically 
in  Fig.  2;  the  branches  belonging  to  negative  p-values  are  excluded.  If  the  wave¬ 
guide  is  lossless  these  complex  conjugate  modes  do  not  contribute  to  the  power 
transfer  in  the  waveguide  (CHORNEy,  1962). 

If  losses  are  introduced  into  the  waveguide  system  the  mode  A  -  B  -  C  -  D 
will  split  up  into  three  different  modes  and  the  branches  A-B,  B  -  C  and  C  -  D 
will  not  be  connected  to  each  other  but  to  the  complex  inodes  so  that  the  mode 
system  will  be  for  instance  A1  -A-B  -  B",  B’  -  B  -  C  -  C"  and  C’  -  C  -  D  -  D* . 
This  means  that  the  plasma  mode  will  split  up  into  a  number  of  different,  sepa¬ 
rated  plasma  modes  depending  on  how  many  branch  points  there  are. 

In  the  following  the  waveguide  system  is  again  taken  to  be  lossless  and  the 
whole  curve  A-B-C-Dis  conveniently  called  the  plasma  mode.  This  mode 
shows  many  interesting  properties .  For  large  8-values  the  mode  asymptotically 
approaches  the  line  =  0  (plasma  resonance)  and  arbitrarily  slow  waves 

exist.  The  plasma  mode  is  of  TE-character  when  p  *  0  and  of  TM-character  when 
p  f* oo.  The  interval  B  -  C  exemplifies  negative  dispersion  i.e.  group  and  phase 
velocity  of  opposite  sign. 

x)  The  explicit  form  of  the  dispersion  equation  to  *  f (p)  is  expanded  in  a 
Taylor  series  around  the  branch  point.  Using  the  relations  *  f (P)Q  and 
[9u)/ap]0  a  O  we  find  that  ho  *  f"(p)Q(Ap)  ,  resulting  in  either  real  or  purely 
imaginary  values  of  (Ap) 
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The  plasma  mode  also  exhibits  a  continuous  transition  from  slow  waves  to 
fast  waves  unless  there  is  a  branch  point  where  the  phase  velocity  equals  the 
velocity  of  light. 

In  Figs.  3  and  4  dispersion  curves  are  shown  for  the  plasma  mode  with  two 
different  surrounding  materials,  viz.  vacuum  (e^  *  1)  and  Pyrex  (er  *  4.5).  In 
these  two  figures  the  diameter  ratio  r^/r 2  is  kept  constant  and  the  dimension¬ 
less  parameter  fu'^r^/c)  is  varied.  To  the  right  in  each  figure  there  is  a  dot- 
dash  curve  representing  the  quasistatic  approximation,  obtained  by  neglecting 
9A/>t  in 

E  =  -  grad  V  -  0A/»t 


For  a  given  p-value  the  exact  solutions  yield  lower  frequencies  than  the 
approximation.  When  {3  -•  <*>  the  exact  and  approximate  solutions  converge.  The 
figures  also  show  that  the  approximation  is  good  only  for  relatively  low  plasma 
densities.  Fig.  5  shows  a  set  of  dispersion  curves  where  u  r_/c  is  kept  constant 

p  i 

and  the  diameter  ratio  has  been  varied.  Together  Figs.  3-5  show  that  the  curves, 
and  especially  the  negative  dispersion,  are  strongly  dependent  upon  the  geometry 
and  the  plasma  density.  The  magnitude  of  the  group  and  phase  velocities  in  the 
region  where  negative  dispersion  exists  is  jv^j  ~  v/IO  and  Jv^J  v/100. 

Using  the  notations  in  Fig.  2,  we  define  the  relative  bandwidth  for  the 
range  of  negative  dispersion 


Aw  _  u>  -  u 
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We  can  see  from  Figs.  3  and  4  that  for  a  given  geometry  the  relative  bandwidth  of 
the  exact  solutions  has  an  upper  limit  given  by  the  relative  bandwidth  obtained 
by  the  quasi-static  approximation.  In  Fig.  6  this  maximum  value  is  plotted  vs. 
the  diameter  ratio;  we  find  that  the  negative  dispersion  disappears  when 
r./r,  >  0.45  and  has  a  maximum  when  r./r  -*  0*  Fig.  6  also  shows  that  a  higher 
gives  a  larger  relative  bandwidth. 


III.  Plasma  column  surrounded  by  a  two-layer  dielectric  enclosed  in  a  metal  tube. 

The  waveguide  system  of  Pig.  1b  consists  of  a  central  plasma  column  surround¬ 
ed  by  two  different,  concentric,  homogeneous  dielectrics  inside  a  metal  tube.  The 
characteristic  equation  can  be  written  in  determinant  form 
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Hie  following  abbreviations  have  been  used: 
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and  similarly  for  F^  and  F^2 
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o 
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Ihe  short  forms  (J  J21  ...)  are  defined  in  the  same  way  as  in  section  II. 

The  principal  behaviour  of  the  dispersion  curves  given  by  eqs.  (5)  is  the 
same  as  for  eq.  (2),  see  Fig.  2.  What  was  said  about  the  plasma  mode  in  section 
II  is  valid  also  for  this  case  except  that  the  position  of  the  asymptote  is  now 

given  by  6  ♦  e  =0.  This  means  that  only  the  inner  dielectric  affects  the 

?  § 

asymptote  and,  as  the  cut-off  frequency  is  less  dependent  upon  the  inner  dielectric, 

a  proper  choice  of  dielectric  material  may  give  rise  to  negative  dispersion. 

Solutions  of  eq.  (5)  are  shown  in  Figs.  7  and  8  for  £  *  4.5  and  c  *  1. 

§  £ 

This  corresponds  to  a  plasma  column  inside  a  glass  tube.  In  Fig.  7  the  dot-dash 

curve  again  represents  the  quasi-static  solution.  Calculations  show  that  as  long  as 

to  r_/c  <  1,  the  approximation  is  in  good  agreement  with  the  exact  solutions, 

P  4 

except  for  small  fi-values .  In  this  case  (Fig.  7)  the  maximum  relative  bandwidth 
(33  %)  is  considerably  larger  than  for  glass  only  (15  %)  or  vacuum  only  (10  %). 
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IV.  Plasma  column  surrounded  by  an  infinite  homogeneous  dielectric. 

The  waveguide  system  of  Fig.  1c  consists  of  a  plasma  coition  in  an  infinite 
homogeneous  dielectric.  The  dispersion  equation  is 


r  *p  VVl2  _  ‘'r  VP2V1  r  1  V^l' 

^h1r1  In®l1r1>  h2r1  Kn(b2T^  ^  ^h1r1  In<h1r1) 


I_ (h-r.) 


1  ^2r  1  ^  1  =  2  2[~  1  _  1  "I5 

h2r1  Kn^h2E1^  ^  ~^hiri^2  (h2t-|)2  J 


(6) 


where  I  ,  K  and  their  derivatives  are  modified  Bessel  functions  of  the  first  and 
n  n 

second  kind  and 


h.  =  k  (t  -  e  ) 

TO  p 

h  =  k  (t2  -  c  j5 
2  o  r 


k  =  ti)/c 
o  . 


t  =  |5/k 


Solutions  to  eq.  (6)  are  Shown  in  Fig.  9  where  the  two  different  families  of 
curves  correspond  to  different  e^.  The  dot-dash  curve  represent  the  quasi-static 
approximation . 

A  comparison  between  Figs.  3,  4  and  7  and  Fig.  9  shows  that  the  surface- 
wave  mode  in  the  waveguide  systems  of  Figs.  1  a-b  is  closely  related  to  the 
surface-wave  mode  on  a  free  plasma  column  -  as  was  stated  in  section  1. 


Experiments 

Measurements  of  the  dispersion  have  been  made,  using  the  experimental  set-up 
shown  in  Fig.  10,  The  waves  could  be  excited  at  either  or  both  ends  of  the  tube. 
Different  methods  of  coupling  microwave  energy  to  the  plasma  mode  have  been 
studied  experimentally  by  KERZAR  (1963).  With  reference  to  his  results  we  chose 
a  single  loop  that  could  be  turned  around  in  the  azimuthal  direction. 

The  plasma  was  the  positive  column  of  a  low-pressure  hot-cathode  mercury 
discharge.  The,  indirectly  heated, oxide  cathodes  were  shielded  from  direct  ion 
bombardment  by  a  metal  screen  and  the  transition  from  the  cathode  region  to  the 
useful  plasma  column  was  made  long  and  smooth.  This  arrangement  gave  very  small 
c«c illations  of  a  low  frequency  (10-100  kHz)  and  the  discharge  seemed  to  be 
stable  in  a  relatively  large  teaperature  interval  (10  °C  -  50  °C).  During  the 
measurements  the  temperature  of  the  coldest  wall  of  the  plasma  tube  was  kept 
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constant  at  20  °C  1  0.5  °C  by  means  of  forced  cooling  and  a  heat  exchanger.  Thus, 

_3 

the  estimated  neutral  gas  pressure  was  1.2*10  mm  Hg.  The  current  through  the 
plasma  tube  was  between  250  -  500  mA  and  stabilized  by  a  regulated  power  supply. 

The  dimensions  of  the  plasma  tube  were  L  =  750  mm,  r^  *  6  mm,  r2  *  7.5  mm  and  of 
the  surrounding  metal  tube  tj-  30  sib  or  16  mm.  Hie  measurements  were  made  at 
frequencies  in  the  range  900-1400  MHz. 

To  measure  the  propagation  constant  of  waves  propagating  in  a  structure  there 
are  basically  two  different  methods  (i)  analyzing  a  standing  wave  pattern  or  (ii) 
analyzing  an  interference  pattern.  The  different  patterns  were  in  this  experi¬ 
ment  plotted  on  an  xy-recorder. 

The  plasma  density  was  measured  by  a  microwave-cavity  method  developed  by 
HEDVALL  (1963),  which  takes  the  surrounding  glass  tube  into  account.  Fig.  11  shows 
a  typical  example  of  the  plasma  density  plotted  against  position  along  the  plasma 
column;  the  density  is  not  constant  along  the  tube.  The  influence  of  this  density 
gradient  upon  the  interference  patterns  was  a  continuous  change  of  wavelength 
(phase  velocity)  along  the  tube.  This  density  inhomogeneity  has  been  observed 
earlier  by  von  ENGEL  (1955),  AGDUR,  KERZAR  and  NYGSREN  (1963),  and  Me  GHOCDf  (1963), 

The  experimental  results  are  shown  in  Fig.  12  together  with  theoretical 
curves  for  (i)  a  plasma  column  surrounded  by  vacuum, (ii)  a  plasma  column 
surrounded  by  Pyrex,  (iii)  a  plasma  column  surrounded  by  a  two-layer  dielectric. 

The  figure  shows  that  the  agreement  between  theory  and  experiments  is  good  when 
the  glass  tube  is  taken  into  account.  The  measurements,  which  were  made  at  diffe¬ 
rent  plasma  densities  (0.5  <  to  r_/c  c  1.5)  and  with  different  geometries  (r,  = 

P  J 

*  30  nan  and  r3  3  16  mm),  were  in  all  cases  in  good  agreement  with  the  theory. 

When  the  phase  velocity  was  equal  to  or  larger  than  the  velocity  of  light  in  the 
surrounding  glass  tube  the  waves  disappeared  and  no  measurements  were  possible. 

At  low  frequencies  another  limitation  prevented  measurements,  viz.  the  symmetrical 
mode  that  was  more  strongly  excited  and  much  less  attenuated  than  the  dipole  mode . 

If  we  consider  the  theoretical  curves  (Fig.  7)  in  connection  with  the  axial 
inhomogeneity  in  the  plasma  density  (Fig.  11)  we  find  that  there  should  exist  a 
cut-off  phenomenon  when  waves  are  propagated  towards  a  lower  plasma  density  (in 
this  case  from  anode  to  cathode).  This  cut-off  was  clearly  visible  in  the 
experiment. 
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The  maximum  relative  bandwidth  as  a  function  of  the 
diameter  ratio.  The  relative  bandwidth  is  defined  in 
Eq.  (4)  and  Fig.  2.  Two  curves  are  drawn  for  different 
values  of  e  . 


frequency  for  ^  &  °f  the  no^alized 

-wfrz.'zr.zr.rr.?'* 

constant.  P  3'  • 2  are  kept 


The  propagation  constant  as  a  function  of  the  normalized  frequency  for  the 
surface-wave  dipole  mode  in  a  waveguide  system  according  to  Fig.  l.c. 
Curves  are  drawn  for  2  different  values  of  with  cu^r ^ / c  as  a  parameter. 
The  dot-dashed  curves  are  the  corresponding  quasi-static  approximations. 
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IV. 


PLASMA  RESONANCE  IN  A  GERMANIUM  ROD 


S3. 


L _ INTRODUCTION 

The  study  of  cyclotron  resonance  in  semiconductors 

has  yielded  much  valuable  information  on  their  energy-band 
12  3 

structure  1  ’  .  In  order  to  increase  resonance  detectability, 

a  high  carrier  concentration  is  desirable,  but  this  introduces 

magneto-plasma  resonance  effects  as  the  carriers  oscillate 
4  5 

collectively  '  .  Earlier,  pure  plasma  resonance  has  been  observed 

in  flat  samples  of  semiconducting  material,  where  it  is  exhibited 

in  a  transition  from  a  highly  reflecting  to  a  transmitted  state,  as 

6 

the  frequency  of  the  incident  wave  is  varied  through  resonance  . 

After  the  investigation  described  in  the  present  paper  was  finished 

(Dec.  1961),  it  became  known  to  the  author  that  magnetoplasma 

resonance  absorption  in  small  Ge  specimens  in  a  constant  field 

had  been  measured  and  successfully  interpreted  theoretically  by 

7 

Michel  and  Rosenblum  .  They  obtained  good  agreement  between 
electron  density  determinations  with  Hall  effect  and  with  magneto¬ 
plasma  resonance  absorption.  In  the  present  paper  the  author 
investigates  plasma  resonance  absorption  with  zero  static  magnetic 
field  for  a  different  geometry,  i.  e.  a  thin,  semiconducting  rod 
situated  in  an  inhomogeneous  field.  This  work  is  connected  with 
similar  studies  on  gas  discharge  plasma®  and  metal  plasma 
carried  out  at  this  laboratory.  The  agreement  between  density 
determinations  with  Hall  effect  and  with  plasma  resonance  absorption 
turns  out  to  be  good  also  in  the  present  case  of  a  rod. 
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2. _ THEORY 

The  theory  presented  here  describes  the  effect  on  the  wave 
propagation  in  a  rectangular  wave-guide  that  is  exerted  by  a  thin  trans¬ 
verse  plasma  rod  perpendicular  to  the  electric  field  vector.  The  wave¬ 
guide  propagates  the  TE^-mode  only  and  the  plasma  is  rather  flossy. 
The  configuration  is  shown  in  Fig.  1. 


2.1. 

Plasma  permittivity 

The  plasma  rod  is  treated  as  having  an  isotropic  permittivity 
(relative  dielectric  constant)  given  by  the  formula^ 

\2 

J  (1) 


(o>  /a>) 

.  =  [e  - 

3  /si* 


(1  +^i-) 


where  e  =  dielectric  constant  of  the  crystal  lattice, 
o)p  =  plasma  frequency, - 

t  =  relaxation  time  of  charge  carriers  in  the  plasma. 

We  consider  an  n-type  semiconductor  of  cubic  symmetry  in 

equilibrium  at  low  temperature.  Thus,  we  may  neglect  contributions 

from  positive-hole  conduction  and  thermal  motion  of  the  electrons.  The 

cubic  symmetry  is  sufficient  to  make  the  permittivity  a  scalar  in  spite 

of  the  anisotropy  of  the  crystal.  Thus,  the  permittivity  formula  will 

hold-  far  n-doped  Ge  and 

2  e  n0 

cu  = - *• 

d  e  -m* 


where  n^  =  density  of  conduction  electrons 

m*  =  effective  mass  of  conduction  electrons. 

There  has  been  some  confu  ion  in  the  past  about  the  so-called 
Lorentz  contribution  to  the  local  electric  field  acting  on  the  electrons1®. 
The  inclusion  of  such  a  contribution  in  the  formula  for  the^purmittivity 
is  discussed  in  the  Appendix. 


2.2. 

Plasma  rod  in  a  waveguide 


The  theory  of  microwave  scattering  by  a  plasma  rod  across 
a  waveguide  is  given  briefly  in  tkesAppandax.  A  thin  rod  of  plasma,  that 
is  lossy,  but  not  too  dense,  is  seen  to  be  adequately  described  by  the 
equivalent  lumped  circuit  in  Fig.  2. 

The  lumped  elements  are  given  by  the  equations 


C  -  1  ~  1 

ci  -r+T 


«0*aro 


£offaro 


At  +  02 


At  + 1  r 


E  o  Tar0 


e0»ar0 


and  the  characteristic  admittance  of  the  waveguide  is  given  by  Eq*  (57). 

The  displacement  current  due  to  vibrating  lattice  electrons 
flows  through  Cj.  The  series  resonance  circuit  (C,  L,  r)  carries  the 
convection  current  of  free  electrons.  Both  these  currents  3hunt  the 
transmission  line.  There  are  thus  two  different  resonances  in  the 
system:  one  series  resonance,  C0j  =  V 1/LC,  giving  a  high  shunt  current 
and  a  reflection  maximum,  and  one  parallel  resonance,  giving  a  low 
shunt  current  and  a  reflection  minimum,  co ^  =V(C  +  Cj)/VLCCj,  where 

®2  >  ®r 

Expressed  in  physical  quantities  one  has  from  Eqs.  (2)  -  (4): 

(«p/»l)2  =(«+!)  (6) 

(“p/^2  =  (*  -  0  (7) 

The  series  resonance,  cu^,  is  called  plasma  resonance  and  is 
pronounced  only  when  the  rod  is  thin  and  the  damping  is  not  too  high. 


The  validity  range  of  the  equivalent  circuit  of  Fig.  2  is  given 
by  four  conditions: 

1.  The  rod  diameter  shall  be  small  compared  to  the  waveguide 
height: 

(2r0/b)2  «  i  (8) 

2.  The  rod  diameter  shall  be  small  compared  to  the  free  space 
wave-length: 

(2jrr0/X)2  «  1  (9) 

3.  The  rod  diameter  shall  be  small  compared  to  the  wave-length 
in  the  rod  (or  the  skin  depth,  if  the  waves  are  evanescent): 

I  *  3  |  <2Tr0A)2«l  (10) 

4.  The  electron  relaxation  time  must  not  be  too  long  in  order 
that  the  dipole  field  shall  predominate  as  is  assumed.  This 
is  expressed  in  the  following  condition: 

(Z'i'oA)2  <  (•■) 


3.  EXPERIMENTAL  SET-UP 


In  order  to  observe  plasma  resonance  in  a  semiconductor  it 
is  imperative  to  reduce  the  damping  by  choosing  the  semiconducting 
material  and  the  experimental  conditions  carefully  with  the  aim  to  get 
an  (<DjT  )- value  as  high  as  possible.  As  cuj  and  are  coupled  by  the 
resonance  condition  (6),  a  high  implies  a  low  resistivity  semicon¬ 
ductor,  which,  in  turn,  would  have  a  small  T  -value.  A  decrease  in 
the  temperature  will  increase  t  ,  but  at  the  same  time  the  electron 


density  and  therefore 


and  (u,  will  decrease. 
P  1 


The  best  choice  avai¬ 


lable  appeared  to  be  a  Ge-rod,  doped  with  Sb  so  as  to  have  a  room- 


temperature  resistivity  of  6.  55  ft  cm,  and  to  put  it  in  an  8  mm  wave¬ 


guide  at  liquid  hydrogen  temperature.  This  should  give  an  (ojjT  )-value 
of  about  3. 


An  analysis,  based  on  the  equivalent  circuit  of  Fig.  2,  shows 
that  the  reflection  from  such  a  rod  is  almost  unaffected  by  the  passage 
through  plasma  resonance.  Therefore  we  chose  to  observe  the  ab¬ 
sorption  resonance,  which  remains  pronounced  down  toaijts  2.  In¬ 
stead  of  varying  the  frequency,  o>,  so  as  to  pass  through  resonance,  we 
varied  the  plasma  frequency  (electron  density),  by  changing  the  rod 
temperature.  The  electron  density  covers  about  three  octaves,  when 
the  temperature  is  increased  from  the  triple  point  to  the  boiling  point 
of  hydrogen  (13.8°K  to  20. 4°K).  This  manner  of  changing  the  density 
of  free  electrons  has  the  advantage  that  thermal  equilibrium  is  conser¬ 
ved.  The  density  is  approximately  the  same  everywhere  in  the  rod. 

The  experimental  set-up  is  shown  in  Fig.  3  and  described 

below; 

The  microwave  power,  46  mW  at  X  =  8.  65  miq  is  fed  from 
a  klystron  through  a  calibrated  attenuator  to  the  section  containing  the 
Ge-rod  (cross-section  0.45  x  0.60  mm).  By  means  of  a  shaft,  the  rod 
can  be  rotated  manually  through  180°  around  its  axis  from  outside  the 
cryostat.  The  line  is  terminated  by  a  matched  load  and  the  waveguide 
section  inside  the  cryostat  is  filled  with  helium  and  it  is  pressure  iso¬ 
lated  from  the  outside  by  means  of  individually  matched  mica  windows. 

The  cryostat  has  an  outer  container  with  liquid  nitrogen  and 
an  inner  container  with  liquid  hydrogen.  The  inner  container  can  be 
evacuated  above  the  hydrogen  surface  to  lower  the  temperature  and 
can  be  filled  with  helium  gas  to  avoid  air  being  sucked  into  the  cryostat. 

The  resistance  of  the  Ge-rod  has  a  large  negative  tempera¬ 
ture  coefficient.  Thus,  the  microwave  power  absorbed  by  the  rod  can 
be  measured  by  connecting  the  rod  to  a  self-balancing  thermistor  bridge 
(Hewlett-Packard  power  meter).  The  rod  is  shunted  with  a  variable 
resistance  so  that  the  bridge  can  be  balanced  for  different  rod  resis¬ 
tances  without  changing  the  built-in  bridge  resistance  of  the  meter. 

When  the  shunt  resistance  is  known,  the  power  absorbed  in  the  rod  may 
easily  be  calculated  from  the  reading  of  the  power  meter. 

The  self-balancing  bridge  thus  keeps  the  rod  resistance 
constant  and  consequently  also  the  density  of  free  electrons.  The 
surrounding  hydrogen  bath  is  kept  at  a  temperature  just  below  that  of 
the  rod.  The  calibrated  attenuator  is  used  to  regulate  the  microwave 
power,  so  that  the  power  meter  shows  a  relatively  constant  deflection, 
in  spite  of  the  variations  in  the  absorption  coefficient  with  temperature. 
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4. 

EXPERIMENTAL  RESULTS 
4,  i .  Measurements 

As  described  above,  the  ratio  of  the  microwave  power  absorbed 
in  the  rod  to  the  power  incident  on  the  rod  -  the  absorption  coefficient  - 
was  measured  as  a  function  of  the  rod  resistance.  The  experimental 
results  are  plotted  in  Fig.  4. 

The  absorption  coefficient  shows  the  expected  resonance  beha¬ 
viour.  The  points  in  Fig.  4  are  measured  with  the  rod  oriented  as  shown 
in  Fig.  5,  which  gives  maximum  absorption  when  the  rod  is  rotated. 
Technological  difficulties  made  it  necessary  to  use  a  rod  with  a  rectan¬ 
gular  cross-section  instead  of  the  circular  one,  on  which  the  theory  is 
based.  Thus,  the  isotropy  of  the  permittivity  could  not  be  checked  but 
the  measured  ratio,  0.  8,  of  the  variation  of  absorption  with  orientation 
of  the  rod  is  not  inconsistent  with  the  assumption  of  isotropy. 


4.  2. Comparison  with  Theory 

The  dimensions  of  the  rod  (length,  1  =  11.2  mm,  cross- section 
2 

S  =  0.  45  x  0.  60  mm  )  determine  the  product  (n^Rp ),  where  R  is  the  dc 
resistance  and  p  the  carrier  drift  mobility  in  the  rod 

)  =  ~e5”  ( 12) 

one  has 

(R  in  «)  (13) 
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p  =. 
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when  to 
P 


An  analysis  of  the  equivalent  circuit  in  Fig. 
is  varied,  the  absorption  maximum  occurs  at 


2  shews  that. 


(14) 


which  should  be  compared  with  the  plasma  resonance  condition  of  Eq.  (6). 
In  the  derivation  of  Eq.  (14)  the  variation  of  t  with  temperature  can  be 
neglected  as  long  as  (cot  )  is  large  enough  to  give  an  observable  absorp¬ 
tion  resonance. 
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Finally  for  t  -  16,  Eqs.  (13)  and  (14)  give 

cot  =  (15) 

»  res 

From  Fig.  4  we  find  the  rod  resistance  at  the  absorption 

maximum  to  be 

R  =  465  a 
res 

Eq.  (15)  now  gives 

(<DT)ree  =  2*  53 

t  s  1.16  •  10~*  *  sec 
res 


There  is  a  slight  inconsistency  in  this  derivation,  as  the 
relaxation  time  in  Eq.  (13)  has  its  dc  value,  while  the  relaxation  time 
in  Eq.  (14)  should  be  taken  at  the  microwave  frequency.  These  two 
values  are  not  necessarily  identical.  As  (cor)2  »  1,  Eq.  (14)  is  almost 
^dependent  of  T  ,  which  means  that  Eq.  (15)  gives  the  dc  value  of  the 
relaxation  time. 

Germanium  has  ellipsoidal  energy  valleys  in  the  conduction 
band,  so  that  *  * 


1 


=  4(_i_  +.L  ,, 

3  vm  u 


(16) 


m  have  been  determined  from 


where  m  ff  ~  1 . 58  •  m  and  mj^  =  0.  082 

cyclotron  resonance  measurements2.  Eq.  (16)  is  valid  if  the  relaxation 

times  along  the  valley  axis  and  perpendicular  to  it  are  identical,  which 
13 

has  been  shown  to  be  a  good  approximation  for  Ge.  Thus,  m*  ~ 

=  0. 12  m  from  Eq.  (16)  giving 

p  =17.0  m2/V 

3.3 


n. 


in!9 
10  m 


sec 

-3 


at  the  absorption  resonance,  which  occurs  at  about  16°K.  These  values 

are  in  good  agreement  with  Hall- effect  determinations  of  the  mobility 

14 

in  similar  Ge  samples 
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To  establish  the  agreement  between  experiment  and  theory 
further,  the  resonance  line  width  will  now  be  investigated.  A  Q- value 
of  the  resonance  is  defined: 


Q  = 

where  index  1  and  2  refer  to  the  points  where  the  absorption  has  half 
its  maximum  value. 


From  the  equivalent  circuit  in  Fig.  2  one  gets 


This  expression  is  valid  only  when  the  absorption  coefficient 
at  resonance  is  much  smaller  than  unity.  It  is  assumed  that 
(orr  =  (an  )reg  (1  +  6)  and  (an  )2  =  (an  )res  (i  -  6), 

where  6  «  1  and  positive. 

2 

If  (an  )  »  1,  Eq.  (l8)is  simplified: 


2(ayr)rea  (yHe. 

[(o>t)2  (o>  /o>)2  -  (a,  (o>  AdJj] 


0)  T 


Q 

~  1  +  Q6 


(19) 


For  our  case,  Eqs  ( 1 3)  and  (17)  yield 


Fig.  4  gives  =  595 £1  and  R2  =  320  £2  so  that 

Q  =  2.98 

If  it  is  assumed  that  6  =  0. 1,  meaning  a  i  10%  variation  of 
t  with  temperature  within  the  resonance  width  (cf.  similar  Ge  specimens 
in  Ref.  14),  Eq.  ( 1 8)  gives  (on  )reg  =  2.22.  This  value,  which  should 
be  valid  at  the  microwave  frequency,  is  about  10%  lower  than  the  dc 
value  computed  from  Eq.  15.  The  accuracy  is  not  good  enough  to 
warrant  any  conclusions,  but  one  should  expect  the  rf  value  to  be  lower 
than  the  dc  value. 

Finally  the  maximum  value  of  the  absorption  coefficient.  A, 
will  be  considered. 


where 


The  equivalent  circuit  in  Fig.  2  yields 
A  A*0* 

max  (1  +  jJ°J/2)Z 

(Q)  (\A)(2»r0/X)2 

A'  . . . -■ .  - - 

(e  +  1)  !  1  -  (X/2a)2  L  / 1  +  (cut)2  -  cut] 


X  =  8.  65  mm  and  a:  7,1  mm  and  b  *  3.  55  mm  are  the  dimensions  of 
the  waveguide. 

From  Fig.  -4-we-obtain  the  value 

A  =  0. 067 
max 

and  this  together  with  the  calculated  value  of  (cut  )  enables  us  to  assign 
an  equivalent  radius  r^  to  the  rod: 

1)  cut  -  2.53  2)  cut  =  2.22 


2rQ  =  0.  76  mm 


2  Tq  =  0.  81  mm 


A  comparison  between  the  equivalent  rod  with  diameter 
0.81  mm  and  the  actual  rod  with  rectangular  cross-section  is  made  in 
Fig.  5. 

It  is  now  possible  to  determine  the  validity  range  of  the 
equivalent  circuit  for  the  Ge-rod  used.  A  plausible  value  of  the  equi¬ 
valent  rod  diameter  is  chosen. 
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2tq  -  0.  65  mm,  cot  =  2.  22  and  (-^E)2  3 

1.  (2r^/b)^  =  0.034  and  Condition  (8)  is  thus  fulfilled. 

2.  (2»r0A)2  =  0.056  and  (9)  "  ”  " 

3.  |e3|  (2rr0/X )  <  1,  if  R  >  250  8.  The  minimum  value  of 
Je3|(2rrQ/\  )2  is  0.37,  which  occurs  atR  =  530  Q.  Cond.  (10) 
is  thus  not  fulfilled  for  R  <  250  Q  and  only  poorly  fulfilled 
for  R  >  25012  (cf.  Fig.  4). 

4.  ^r(to  t  )  ~  =  9  and  Condition  (11)  is  thus  fulfilled. 

In  view  of  the  fact  that  Condition  (10)  is  poorly  fulfilled  for 

our  rod,  the  agreement  with  the  simple  theory  is  surprisingly  good. 

The  determination  of  (cur  )  from  R  ,  using  Eq.  (15),  or 

ires  ires 

from  Q,  using  Eq.  (18),  is  not  influenced  by  an  inclusion  of  a  Lorentz 

contribution  to  the  local  electric  field  acting  on  the  free  carriers  (see  the 

App.,  especially  Eq.  (24)).  The  drift  mobility  would  be  given  by 
(  2  ^  c  )  eT 

p  -  ■  3 — -  and  take  a  value  6  times  larger  than  before,  and  the 

case  is  exactly  the  same  for  the  drift  mobility  determined  from  a  Hall 
effect  measurement.  Those  two  measurements  together  are  thus  not 
sufficient  to  settle  the  question,  whether  a  Lorentz  term  should  be 
included  or  not.  As  mentioned  in  the  App.  a  pulse  drift  measurement 
is  needed  combined  with  either  a  plasma  resonance  or  a  Hall  effect 
measurement.  Such  a  comparison  with  the  directly  determined  drift 
mobility*  ®  shows  that  the  Lorentz  contribution  is  negligible. 


V 


Thera  is  no  indication  of  discrepancy  between  the  simple 
theory  and  the  observed  width  and-shape  of  the  absorption  resonance, 

8 

as  in  the  case  of  polarized  resonant  emission  from  a  mercury  discharge  . 
This  may  be  due  to  the  fact  that  the  theoretical  idealization  of  uniform 
plasma  density  and  negligible  temperature  velocity  is  much  closer  to 
the  truth  for  our  Ge  rod  than  for  a  gas  discharge  in  a  glass  tubing. 

Substitution  of  (cot)  from  Eq.  (15)  into  Eq.  (13)  shows  that 

2 

ojp  and  therefore  also  the  electron  density  is  directly  determined  by  the 
rod  resistance  at  resonance.  As  already  observed  by  Michel  and  Ro- 
senblum^,  the  plasma  resonance  absorption  thus  could  be  used  to  de¬ 
termine  the  charge  carrier  density  in  semiconductors  directly  without 
recurrence  to  Hall  effect  measurements  with  their  uncertainty  a3  to 
the  theoretical  value  of  the  ratio  drift  mobility  to  Hall  mobility.  The 
application  and  accuracy  of  the  method  are  limited  by  the  low  Q- values 
attainable. 


APPENDIX 


(TiC  looale^cctric  field 

There  is  a  theoretical  possibility  that  the  local  electric  field, 
E^qc<  acting  on  the  conducting  electrons,  is  different  from  the  macro¬ 
scopic  field,  E,  in  the  crystal.  In  this  general  case  the  permittivity 
for  the  crystal  with  plasma  is 


=  [e  -  (- 


(“  A>) 


(1  +-pi-  ) 

\  1CJDT  ' 


At  the  lattice  points  in  a  crystal  of  cubic  symmetry  there  is  a 
local  field  ^ 

E,  =-^-E  (24) 


This  field  could  possibly  be  the  effective  field  for  the  conduc¬ 
tion  electrons  also.  There  are,  however,  strong  reasons  to  discard 
this  assumption.  It  would  mean  that  the  theoretical  value  of  the  ratio 
drift  mobility  to  Hall  mobility  should  be  multiplied  by  the  facto 
in  contradiction  to  the  general  agreement  between  measured  ratio 

values  for  Ge  and  Si  and  the  theoretical  values  based  on  E.  =  E.  The 

loc  if 

drift  mobility  can  be  measured  directly  by  the  pulse  drift  method  and 
the  Hall  mobility  by  a  Hall  affect  measurement.  It  should  be  observed 
that  ^  ^  E  =  6  for  Ge. 

It  is  also  possible  that  the  local  field  is  influenced  by  the  dipole 

moment  of  the  vibrating  conduction  electrons  themselves,  as  has  been 
12 

proposed  .  In  considering  this  possibility  we  do  not  include  die  effect  of 
a  contribution  from  that  source  to  the  local  field  at  the  crystal  lattice 
points.  We  get 


E.  =  [  1  + 
loc 


1/3  (m  /co) 


r1. 


(ffl  /cd) 

e  3  =  {  e  -  - , - E__ - j — -  }  (26) 

li+±(*pM *  +  -±7] 

Thus  Re(c  3)  0  could  never  be  reached  for  an  s  larger  than  3. 

This  is  contradicted  by  the  observed  rapid  change  from  transmission  to 
total  reflection  in  heavily  doped  germanium  in  die  far  infrared  ,  which 
is  believed  to  be  a  consequence  of  passage  through  Re^j)  -  0. 
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Theorfr^f  h  ttircular  plasma  rod  across  a  waveguide 

The  emact  treatment  of  a  circular  plasma  rod  across  a  wave¬ 
guide  is  indicated  below  (see  Fig.  1),  but  only  a  limiting  case  valid  for 
a  thin  rod  with  comparatively  short  electron  relaxation  time  is  worked 
out  completely.  The  fields  have  the  time  dependence  exp(imr  )  and  the 
temperature  is  so  low,  that  the  temperature- induced  electrostatic  field 
solution  in  the  rod  may  be  neglected. 

Only  the  field  vectors  in  the  axial  direction  (r- direction)  are 
given,  as  the  transverse  components  are  easily  deduced  from  the 
formulas: 

-  BE  A 

f> *,  *  Vur*  +  k^o<z  *  7tHz>  <i7> 

=  (j)2J(£*VtEz)4  (28) 

where  p  e  ^|(S)2  -  (£)2 

z  is  die  unit  vector  in  the  z-direction.  Index  t  stands  for  "transverse11 
to  the  z-direction.  Formulas  (27)  and  (28)  are  applicable  outside  the 
rod;  inside  the  rod  p  is  replaced  by  k,  defined  as 


The  TEjg-mode,  incident  in  the  y-direction  and  arbitrarily 
normalized,  has  the  form: 


Ejp  =  0  (29) 

H<l)»  cos  (£)  exp  (-ipy)  (30) 


Using  the  expansion  exp(-i0y)  =  2  J  (dr)  exp(-in<p),  where 

n  n 

y  =  r  •  sintp  and  the  summation  is  from  -  oo to  +  co,  changes  this  into 
E^  =  0  (31) 

H«  =  g  cos(^)  Jn(Pr)  exp(-in(?)  (32) 

Although  E^  =  0,  the  electric  field  inside  the  rod  and  the  scattered 
electric  field  will  have  E  /  0.  This  field  component  is  excited  by  the 
longitudinal  current  in  the  plasma  rod  induced  by  the  H^- component. 
Inside  the  rod  we  get 


4 


(33) 

(34) 
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■  g  »in(^)  AnJn(kr)  exp(-in<J>)  (33) 

*  n  co*(?)  cnJn(kr)  exp(-inq>)  (34) 

If  the  walls  at  x  *  -  b/2  were  moved  out  to  infinity,  the 
scattering  fields  outside  the  rod  would  be  represented  by  simple  cylind¬ 


rical  waves 


=  2  sin(££)  Bn  H^(Pr)  exp(-inf) 
=  2  cos(~)  D  H^(Pr)exp(-in<p) 


z  n  *  a  '  n  n  *  '  ’  '  ’ 

The  walls  will  introduce  mirror  images  of  the  rod  and  images 
of  images  in  an  infinite  array.  Thus,  die  boundary  conditions  at  die 
walls  are  fulfilled  by  the  complete  scattering  fields: 

sin(-)B  H^(pr  )  exp(-in<p  )  (37) 

z  nm  'a,nnvrm/rv,m/  '  ' 

=  22  cos(— )  D  (Pr  )  exp(-in<p  )  (38) 

z  nm  '  a  '  n  n  m'  m7  '  ' 

where  r  and  <p  are  defined  in  Fig.  6.  It  should  be  observed  that 
mm 


A  =  -A 
-n  n 

C  =  G 
-n  n 


B  *  -B 
-n  n 

D  =  D 
-n  n 


For  a  point  in  the  waveguide  far  away  from  the  rod,  |  y  |  »  b, 
it  is  possible  to  deduce  the  following  asymptotic  expressions: 

E^  -  0  (39) 

H  zZ)  -  g  VG  cos(lT)  (7yf)n  Dn  exP(_ip  M)  (40) 

which  are  valid  as  long  as  the  waveguide  can  propagate  only  the  TE^q- 
mode.  The  scattered  wave  may  thus  be  described  as  one  back-scatte¬ 
red  wave  giving  a  voltage  reflection  coefficient 

r=-„|f  <-‘>X  HD 

and  one  forward- scattered  wave,  which  added  to  the  incident  wave 
gives  a  voltage  transmission  coefficient 

T*<ltin;  £D„>  <«> 

r  and  T  are  referred  to  y  =  0  (a  plane  through  die  centre  of  the  rod), 
and  give  a  complete  description  of  the  effect  of  the  rod  on  die  trans¬ 
mission  line.  The  D  can  be  determined  from  the  conventional  boun- 

n 

dary  conditions  for  die  fields  at  the  surface  of  the  dielectric  rod. 
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(pr0)2«  landlkr0JZ«l 

The  case  Otp)2  «  1  and  |kr0  f  «  1  will  new  be  considered 
in  detail.  As  a  zeroth  order  approximation,  the  contribution  from  the 
image  fields  to  the  fields  at  the  rod  boundary  is  neglected.  We  get  to 
the  lowest  order  in  (0rQ): 

D0>  =-‘7<&>Z£2>4<‘3-,>  <«> 

d<o)  -  -  mi  &2n  *«  •  *  -  »..■».  <«> 

The  values  of  are  true  at  Re(c^)  »  - 1  only  if 

f§  £^>2  <«> 

For  a  semiconductor  this  is  equivalent  to  the  condition 

(ojt)  «  PKfJ'-1.)  (46) 

2W-/)2 

Only  Dj  and  D  ^  will  contribute  to  r  and  T  to  the  lowest 
order  in  (PTq): 

0>r„  ,  t  ,  -  1 

(47) 


r(°)  £  (”0.)2  i3'  N 

=  i  + 


(48) 


The  restriction  on  Im(c  in  Eq.  (45)  is  equivalent  to  stating 
that  Jr^|  «1  is  a  condition  for  T  ss  to  be  valid.  A  further  step 
in  the  approximation  procedure  will  replace  this  condition  by  the  less 
restrictive 


u>r 


0X2 


(49) 

(50) 


fMe3)|  >  7  {——) 

which  gives:  (tor)  ,<  — 

[o^q/c) 

The  new  condition  guarantees  that  still  only  the  dipole  mode, 
Dj  and  D  j,  will  contribute  to  the  reflection  coefficient  at  plasma 
resonance,  Re(c  j)  =  -1,  If  |lm(c^)j  is  further  decreased,  at  first  the 
quadrupole  mode,  and  D  and  later  the  sextupole  mode,  and 
D_j,  and  higher  modes  contribute  to  the  absorption  and  reflection  at 
plasma  resonance. 
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In  this  first  order  approximation,  the  image  fields  are  added 
to  the  fields  at  the  rod  boundary  and  to  this  order  of  approximation: 

D(0) 

DW  *  - . . r-*“iwr  (51) 


“  '  *  D'i 


which  gives 

r(0s 


[i  -rW) 

T(,)=1,tr<1))  -  t,  .  *rW) 


The  result  is  thus 


“  “u"r(o) 
r”  tTTTWT] 
t*[i  +  r]  M  - 


U  -  r'0'] 


valid  for  (0rn)  «  1 


with  the  addition  of  condition  (50). 

The  form  T  «  (l  +  T)  tells  us,  that  the  rod  constitutes  a  lumped 
shunt  element  in  the  waveguide.  It  is  easily  observed  from  Eq.  (54), 
that  the  normalized  shunt  admittance  is  given  by 

^=-2r(0)  (56) 

c 

where  the  characteristic  admittance  of  the  waveguide,  Yc,  is  given  by 
Yc 

Using  Eqs.  (47),  (l)  and  (57)  the  shunt  admittance  of  Eq.  (56) 


caui  be  written 


Y  =  {i 


TT}  + 


(■  +  *r +  i<D  l«  .t  ‘X  trL  ii-Lil 

7T  ,  2  u>  2 


-]  (58) 


2cd  t 
P 


A  suitable  equivalent  circuit  description  of  the  rod  in  the 
waveguide  is  now  obvious  and  given  in  Fig.  2  with  die  lumped  elements 
a^  in  Eq.  (2)  through  (5). 
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Thin  plasma  rod  mounted  across  a  rectangular 
waveguide  parallel  to  the  broad  walls. 


Equivalent  lumped  circuit  describing  a  thm  plasma 
rod  across  a  waveguide. 
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Fig.  3;  Block  diagram  of  experimental  setup  for  plasma 

resonance  measurement  in  Ge. 


5.  Orientation  in  the  microwave 
field  of  actual  rectangular 
rod  and  for  comparison 
an  "equivalent"  circular 

rod.  I 

I 

Power  in  / 


6.  Plasma  rod  in  a  waveguide 
and  its  reflected  images 
in  the  waveguide  walls. 
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V.  MAGNETO- PLASMA  RESONANCE 
Introduction 

If  a  metal-  or  semiconductor-plasma  sample  is  immersed  in  a  d.c.  magnetic 
field,  and  irradiated  by  an  electromagnetic  wave  with  a  frequency  much  less 
than  the  plasma  frequency  and  the  gyrofrequency,  then  the  plasma  sample  will 
be  resonant  for  certain  values  of  the  frequency.  These  resonances,  which  are 
commonly  called  magneto-plasma  resonances,  have  been  observed  experimentally 
in  both  metals^’ 2^  and  semiconductors3^.  At  our  laboratory  we  have  studied 
the  associated  propagating  waves  in  a  long  semiconductor  rod. 

The  resonances  observed  can  generally  be  divided  into  two  groups: 

1)  the  main  resonances  with  different  mode  number  in  the  direction  parallel 
to  the  magnetic  field,  and  2)  the  secondary  Cor  satellite)  resonances, 
corresponding  to  different  transverse  mode  number. 

The  main  resonances  can  be  fairly  well  explained  by  assuming  plane  waves 
propagating  along  the  magnetic  field.  (This  simple  approach  is  valid  if  the 
transverse  dimensions  are  much  larger  than  the  wavelength  along  the  magnetic 
field  lines.) 

The  secondary  resonances  (which  yet  have  been  observed  only  in  metal 
samples)  can  not  be  explained  by  the  plane-wave  approximation.  To  obtain  a 
theory  for  these,  one  has  to  solve  the  proper  boundary-value  problem,  lhis  is 
done  in  an  approximate  manner  for  circularly  cylindrical  samples  in  the  theory 
below. 


Theory 

The  ordinary  theory  for  a  magnetized  plasma,  describes  the  plasma  as  an 
anisotropic  dielectric  medium,  having  the  tensor  dielectric  constant: 
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and  er  *  static  dielectric  constant 

u>  ■  signal  frequency,  co  =  plasma  frequency 

P 

uc  *  gyro  frequency  and  v  =  collision  frequency 

provided  the  d.c.  magnetic  field  coincide  with  the  z-direction. 

The  solutions  of  Maxwell's  equations  for  plane  waves  give  two  possible 
normalized  propagation  constants 

^  '  (C/V>2  *  ‘l  1  e2  <2) 


where  v  h  =  the  phase  velocity 
and  c  =  velocity  of  light 

If  now  a)  '»■  |u)  -  jv|  and  >  £  oa  (as  is  generally  the  case  in  the 
c  1  1  p  r  c 

magnetoplasma  resonance  experiments)  only  the  left  hand  circularly  polarized 
wave  Can  propagate  and  the  normalized  propagation  constant  can  in  this  case 
approximately  be  written: 


e  ♦  — ■— 

r  ui*u> 


For  a  metal  plasma  tir/uxu^  and  hence  the  main  resonant  frequencies 
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where 


z  =  length  of  the  sample  in  the  d.c.  magnetic  field  direction, 
p  =  axial  mode  number. 


We  will  now  consider  a  circularly  cylindrical  plasma  sample,  with  radius  a 
and  axial  length  z,  and  with  the  d.c.  magnetic  field  parallel  to  the  axis. 

According  to  AGDUR*^  the  z-  and  <f— components  of  the  fields  inside  the  plasma 
(r  <  a)  can  be  written  (omitting  the  factor  expfjfwt  -  |3z  ♦  n?)]  ): 
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Jn  is  the  n:th  order  Bessel  function  and  is  its  derivative;  A  and 
B  are  amplitude  coefficients. 

Ihese  fields  must  satisfy  certain  boundary  conditions  at  the  plasma  boundary 


r  *  a,  which  gives  the  dispersion  relation.  We  will  here  distinguish  between 
two  cases:  A)  the  plasma  surrounded  by  air,  which  applies  to  most  experiments 
with  metal  plasmas,  and  B)  the  plasma  surrounded  by  a  metal  tube  (d 
which  applies  to  most  experiments  with  semiconductor  plasmas.  These  two  cases, 

0 

which  will  give  different  boundary  conditions,  are  treated  separately  below. 
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Case  A:  Plasma  surrounded  by  air 

In  this  case  the  fields  outside  the  plasma  (r  >  a)  are  assumed  to  be  the 
ordinary  slow  wave  vacuum  fields  in  cylindrical  coordinates,  i.e.: 

ez  =  acn(pc) 

-jcB  3  EK  Cpr ) 
z  n 

n  1  < 6> 

E  3  -  CK  (Sr)  *  ~  EK'  (pr) 

9  pc  n  y  t  n  r 

-jcB  =  -  CK'  (pr)  -  7T  EK  (Sr) 
tp  t  n  H  pr  n  H 

where  p  3  -rpo 

K  and  K1  are  the  modified  second  kind  Bessel  function  of  order  n  and  its 
n  n 

derivative.  C  and  E  are  amplitude  coefficients. 

Matching  the  fields  inside  and  outside  the  plasma  at  the  boundary  r  3  a 
requires  that  the  following  determinant  is  zero: 
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and  if  further  t  is  assumed  to  be  near  the  plane  wave  solution  (i.e.  t  of  the 
order  of  j  e  j )  then: 
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With  these  approximations  one  can  show  that  an  approximate  solution  to  Eq. 
(7)  is 


or 


J„  (m,a)  ~  0 
n  2 


m.a  x 
2  —  nm 


(9) 


where  x  is  the  m:th  rooth  of  J  (x)  =0  (n  ^  0). 
nm  n 

To  obtain  a  simple  expression  for  the  resonant  frequencies  of  a  plasma  rod 

2 

with  the  length  z,  we  approximate  m„  further: 
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This  gives  the  resonant  frequencies: 
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This  equation  gives  both  the  main  resonances  (n  =  m=1,  p=1,  2....) 
and  the  secondary  ones  (p  fixed  and  n,  m  »  1,  2 

If  (x^  z)/(pxa)  1  then  the  plane  wave  approximation  (Eq.  (4))  gives  the 
main  resonances  with  good  accuracy,  and  the  secondary  resonances  will  be  very 
close  to  the  main  ones.  The  secondary  resonances  will  only  be  detectable  if  the 
Q-values  of  the  resonances  are  of  the  order  of  fQ/2Af  or  larger,  where  fQ  is  the 
main  resonant  frequency  and  Af  is  the  frequency  spacing  between  the  main  and  the 
first  secondary  resonance. 

Case  8:  Plasma  surrounded  by  a  metal  tube. 

In  this  case  the  fields  must  satisfy  the  boundary  conditions: 


£  =  E  =  0  for  r  =  a 

z  <P 

With  the  expressions  for  and  E  in  Eq.  (5)  these  conditions  lead  to  the 
following  characteristic  equation  determining  the  propagation  constant: 

Xl“2  -  %  '  Pi  -  “l  (,1> 

where  and  are  given  in  Eq.  (7). 

If  we  introduce  the  same  approximation  as  in  case  A,  Eq.  (11)  can  approximate¬ 
ly  be  written: 
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The  solution  to  Eq.  (11a)  is: 


WB2a)  =  ° 
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where  x  _  is  the  m:th  rooth  of  J  „(x)  =  O;  n  /  0.  From  this  approximate 
n+-1,m  n+1 

equation  one  can  derive  the  propagation  constant  t: 
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Experiments 

We  have  at  our  laboratory  made  an  experiment  with  a  semiconductor  plasma 
in  a  strong  magnetic  field.  The  semiconductor  (at  liquid  nitrogen  temperature) 
was  immersed  in  a  circularly  cylindrical  waveguide,  the  axis  of  which  was  along 
the  magnetic  field  lines.  The  purpose  of  this  experiment  was  to  measure  the 
propagation  and  the  attenuation  constants  of  a  circularly  polarized  wave  in  the 
semiconductor,  as  a  function  of  the  magnetic  field  strength.  The  experimental 
conditions  were  such  that  the  plane-wave  approximation  could  be  used.  In  this 
case  the  ordinary  magneto-plasma  theory  gives  the  propagation  constant  (f3)  and 
the  attenuation  constant  (re)  of  the  extraordinary  wave  (see  Eq.  (2)): 


These  formulas  are  valid  if  wc  »  |u  -  jv| .  The  ordinary  wave  is  strongly 
attenuated  and  will  not  be  considered  here. 

21  3 

The  semiconductor  used  was  indium  antimonide  with  N  ~  10  electrons/m  , 

\x  Cl  30  m  /Vs  and  e  “16.  With  a  signalfrequency  of  35  GHz  and  a  magnetic  field 
strength  of  1  Wb/m^  (*  104  Gauss),  the  above  formulas  give  p  *  |30*10  which 
corrensponds  to  a  wavelength  in  the  staple  of  0.87  am.  The  corresponding  value 
of  a  is  104  m“\  i.e.  an  attenuation  length  of  about  1  cm. 
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The  experimental  setup  is  shown  in  Fig.  (1).  The  incident,  linearly 
polarized  TE^-mode  in  the  empty  waveguide  can  be  divided  into  one  ordinary  and 
one  extraordinary  circularly  polarized  wave.  Only  the  extraordinary  wave  can 
propagate  through  the  semiconductor  and  reach  the  magic  T.  This  signal  interferes 
with  a  signal  passing  through  the  attenuator  and  the  directional  coupler. 

The  magnetic  field  is  obtained  by  discharging  a  capacitor  through  the 
solenoid  around  the  sample.  The  field  reaches  its  maximum  value  (~  1.3*10  Gauss) 
in  0.005  sec. 

An  oscilloscope  trace  of  the  detector  output  as  a  function  of  the  magnetic 
field  strength,  is  shown  in  Fig.  (2).  The  peaks  of  the  curve  are  related  to  the 
wavelength  in  the  sanple,  but  it  is  not  possible  to  get  an  absolute  value  of  the 
wavelength;  only  the  relative  change  of  the  wavelength  can  be  measured.  However, 
by  taking  one  of  the  measured  points  as  reference  point,  one  can  compare  the 
measurements  with  the  theory.  In  a  similar  way  the  measured  attenuation  can  be 
compared  with  the  theoretical  value.  (Fig.  (3).) 

In  our  experiment  the  length  of  the  sample  was  several  (10  -  15)  wave¬ 
lengths  long.  Consequently,  owing  to  the  relatively  strong  attenuation,  the 

3) 

magneto  plasma  resonance  observed  by  LIBSCHABER  and  VEILEX  could  not  be 
observed  in  our  case. 

We  are  presently  preparing  a  new  experiment  in  which  we  hope  to  see  both  the 
main  and  the  satellite  magneto  plasma  resonances  in  indium  antimonide. 

For  dipole  modes  Eq.  (12)  gives 
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m  •  1,  for  the  main  resonance  and 
m  *  2,  3  ...  for  the  satellite  resonances. 

A  numerical  analysis  of  a  reasonable  case,  shows  that  the  spacing  between 
the  main  and  the  proper  secondary  resonances  is  very  small,  but  that  the  un¬ 
loaded  Q-values  of  the  resonances  are  so  large  that  the  different  resonances 
shall  be  separable.  This  means  that  the  coupling  to  the  semiconductor  sample 
must  be  very  small,  in  order  not  to  lower  the  Q-values  too  much.  We  hope  to 
achieve  this  low  coupling,  simply  by  applying  metal  irises  to  the  end-faces  of 
the  sample. 
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VI.  VARIATION  OP  PLASMA  PARAMETERS  ALONG  TOE  POSITIVE  COLUKW  OP  A  MERCURY 
DISCHARGE 

This  investigation  is  a  continuation  of  the  work  described  in  ref.  1,  where 
experimental  studies  of  the  variation  of  the  electron  density  along  the  positive 
column  of  a  low-pressure  mercury  discharge  are  described.  Ihe  electron  density 
was  measured  by  means  of  the  microwave  cavity  method.  It  is  the  purpose  of  this 
experiment  to  obtain  more  information  on  the  axial  variation  of  the  plasma  para¬ 
meters  by  introducing  Langmuir  probes,  so  that  also  the  electron  temperature  and 
the  axial  d.c.  electric  field  in  the  plasma  can  be  measured. 

The  discharge  tube  (Fig.  1)  has  a  length  of  1  m,  an  inner  diameter  of  12  mm 
and  an  outer  diameter  of  15  mm.  The  neutral  gas  density  is  controlled  by  keeping 

the  glass  wall  around  the  cathode  at  constant  temperature  T  ,  which  is  below  the 

S 

temperature  of  the  rest  of  the  tube.  Seven  Langmuir  probes  are  inserted  along 

the  tube  with  a  distance  of  14  cm  between  them.  Only  the  circular  end  surfaces 

are  exposed  to  the  plasma  and  are  situated  on  the  axis  of  the  tube.  A  sawtooth 

voltage  with  a  sweep  frequency  of  50  c/s,  is  applied  to  the  probes  and  a 

logarithmic  current  vs.  linear  voltage  characteristic  displayed  on  an 

oscilloscope  (Fig.  2)  after  compensation  for  the  ion-current.  The  circuitry  is 

2) 

essentially  the  same  as  described  by  HARP  .  The  linear  slope  of  the  characteristic 
is  a  measure  of  the  electron  temperature  and  is  determined  by  comparison  with  a 
reference  line,  the  slope  of  which  can  be  varied.  The  electron  density  can  be 
determined  from  the  saturation  current  when  the  area  of  the  probe  and  the 
electron  temperature  are  known. 

Figs.  3  and  4  show  the  electron  temperature  and  density  at  one  place  in  the 
tube  as  a  function  of  T^  for  various  tube  currents.  The  density  is  measured  with 
the  cavity  method.  The  results  shown  in  the  figures  are  in  qualitative  agreement 
with  KLAP.FELD3? 

The  variation  of  electron  temperature  along  the  tube  is  shown  in  Fig.  5. 

Apart  from  the  maximum  in  the  centre  of  the  tube  the  variation  indicates  that 

there  is  a  small  positive  pressure  gradient  towards  the  anode.  Such  pressure 

4) 

gradients  have  been  discussed  by  EXARFELD  and  POLBTAEV  .  The  maximum  of  electron 
temperature  is  however  not  understood. 

The  determination  of  electron  density  with  probe  measurements  depends  on  the 
effective  area  of  the  probe  which  is  difficult  to  determine  accurately.  The 
relative  areas  of  the  probes  are  found  by  measuring  the  tensity  with  the  micro¬ 
wave  cavity  method  at  a  small  discharge  current  (100  mA).  While  the  probe 
measurements  give  the  density  at  the  axis  of  the  tube,  the  cavity  method  gives 
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approximately  the  mean  density  across  the  tube.  If  we  assume  that  the  radial 
density  distribution  is  constant  along  the  tube  for  I  *  100  mA,  the  relative 
areas  can  be  determined  with  good  accuracy.  Furthermore,  we  assume  that  the 
effective  areas  are  independent  of  the  tube  current  in  the  current  region  of 
our  measurements. 

Fig.  6  shows  the  axial  variation  of  electron  density  measured  by  the  probe 

method  and  the  cavity  method.  The  probe  measurements  are  relative  in  the  sense 

that  N (probe)  is  put  equal  to  N(cavity)  for  a  tube  current  of  100  mA.  The 

difference  between  the  results  of  the  probe  measurements  and  the  results  of  the 

cavity  measurements  indicates  that  the  radial  density  distribution  changes 

along  the  tube.  This  effect  is  more  clearly  demonstrated  in  Fig.  7  where  a  = 

=  N(probe)/N(cavity)  is  plotted.  This  ratio  which  is  normalized  to  unity  at 

100  mA  is  approximately  proportional  to  N  /N  under  the  above  mentioned 

max  mean 

assumptions.  It  is  seen  that  a  is  approximately  constant  over  the  greater  part 
of  the  tube  but  decreases  near  the  anode  and  near  the  cathode.  In  general  it  is 
increasing  with  the  current. 

If  we  assume  that  the  pressure  gradient  causes  the  variation  of  electron 
density  and  temperature,  and  that  these  are  determined  by  the  local  pressure, 
we  can  determine  the  pressure  gradient  from  the  measurements  of  the  electron 
density  distribution  or  the  electron  temperature  distribution  as  plotted  in 
Fig.  8.  If  we  neglect  the  maximum  of  electron  temperature,  the  variation  of  both 
the  density  and  the  temperature  is  in  agreement  with  a  positive  pressure  gradient, 
which  can  be  determined  to  be  about  1.3  x  io-5  mm  Hg/cm.  This  figure  is  in  good 
agreement  with  Klarfeld  and  Poletaev,  who  do  not,  however,  take  into  account  the 
gradient  in  electron  density  in  their  theoretical  considerations. 

The  axial  electric  field  has  also  been  measured  and  was  determined  by  using 
the  probe  characteristics  to  find  the  plasma  potential  at  the  different  probe 
positions.  The  electric  field  was  found  to  increase  slightly  towards  the  anode. 
Similarly,  a  higher  gas  pressure  increased  the  field  at  a  fixed  position  (Fig.  9). 
No  dependence  on  the  current  could  be  observed  (current  range  100  -  600  mA). 

As  a  conclusion  we  can  say  that  the  plasma  parameters  vary  along  the  tube 
approximately  in  accordance  with  an  increasing  pressure  towards  the  anode,  with 
one  major  exception  in  the  electron  temperature.  We  intend  to  extend  the 
measurements  to  higher  gas  pressures  and  possibly  measure  the  pressure  and  its 
variation  along  the  tube  by  spectroscopic  means  in  order  to  get  a  unified  picture 
of  the  phenomena . 
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for  various  tube  currents. 
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Fig.  9, 


Variation  of  the  potential  along  the  tube  for  two  different 
gas  pressures. 
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VIL  EX PERIMELMTS-  OMNOISE  RADIATION- FROM-' A 

PLASMA  IN  A  WAVEGUIDE 

Introduction 

The  work  reported  here  is  a  continuation  of  our  earlier 
experiments  on  noise  radiation  and  coherent  scattering  from  a 
low-pressure  plasma  column  in  free  space  and  in  a  waveguide. 

The  experiments  are  connected  to  the  theoretical  investigation, 
which  is  described  elsewhere  in  this  report.  The  experimental 
investigation  is  being  continued  and  the  results  reported  here 
are  only  preliminary.  AGDUR  et.  al.  have  published  experimental 
data  on  noise  radiation  and  scattering  from  a  plasma  column  in 
free  space  and  in  a  Waveguide.  Strong  polarization  of  noise  was 
observed  in  the  region  of  secondary  peaks.  The  aim  of  this 
report  is  to  present  more  experimental  data  on  noise  radiation 
from  a  cylindrical  plasma  column  situated  in  a  waveguide. 


The  noise  power  radiated  from  the  plasma  column  was 
measured  by  means  of  a  radiometer,  operating  at  a  frequency  of 
3100  Mc/s  with  a  bandwidth  of  10  Mc/s.  A  block  diagram  of  the 
radiometer  used  is  shown  in  Fig.  1.  The  noise  radiation  from 
the  plasma  column  was  compared  with  the  noise  radiation  from  a 
calibrated  noise  source.  In  this  way,  drift  and  nonlinear  effects 
of  the  post  r.f.  amplifier  were  eliminated.  The  combined  noise 
signal  was  then  amplified  in  a  low  noise  travelling  wave  tube 
and  an  I.  F.  amplifier  operating  at  30  Mc/s  with  a  bandwidth  of 
10  Mc/s.  The  amplified  signal  was  rectified  by  a  "square  law" 
detector  and  the  resulting  low  frequency  signal  was  fed  into  a 
selective  amplifier  (f  =  70  c/s,  Af  =  5  c/s)  followed  by  a  synchronous 
detector  and  integrator  with  an  integration  constant  of  15  sec.  The 
lowest  detectible  temperature  change  was  about  20°K. 

The  noise  power  radiated  from  a  plasma  tube  in  a  frequency 
interval  Af  is  spread  over  an  infinite  number  of  modes.  Only  the  noise 
power  carried  by  the  dominant  mode  of  the  waveguide  will  reach  the 
radiometer.  It  should  be  pointed  out  that  no  attempt  was  made  to 
match  the  plasma  tube  to  the  waveguide  system.  Therefore  one 
cannot  relate  the  electron  temperature  and  the  noise  power  in  a 
simple  way.  Simultaneous  measurements  of  the  frequency  spectra 
for  noise,  scattering  and  absorption  were  made  for  two  typical 
positions  of  the  plasma  tube:  one  with  the  E-vector  of  the  dominant 
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waveguide  mode  (TEqj)  parallel  with  the  plasma  tube  and  the 
other  with  the  E-vector  perpendicular  to  the  plasma  tube.  We 
will  here  call  them  transverse  and  longitudinal  polarization. 

See  Fig.  3. 

Two  experimental  mercury  vapour  discharge  tubes  with 
hot  cathodes  were  used,  having  a  diameter  of  7  and  12.5  mm 
respectively.  The  neutral  gas  density  could  be  varied  from 
5  •  lO1^  to  3'  1 020  molekules/m^.  The  electron  density  was 
measured  by  means  of  the  microwave  cavity  method. 

Measurements 

By  means  of  the  experimental  arrangement  the  following 
data  can  be  obtained  as  functions  of  the  electron-  and  neutral 
gas  density:  (i)  radiated  noise  power  from  the  plasma  tube,  (ii) 
absorption,  transmission  and  reflection  of  a  coherent  microwave 
signal  (iii)  exitation  and  propagation  of  slow  wave  plasma  modes 
along  the  plasma  column. 

Transverse  polarization 

The  measured  emission  and  absorption  spectra  have 
substantially  the  same  shape.  This  similarity  indicates  that 
the  plasma  investigated  is  in  thermal  equilibrium. 

The  spectra  of  noise  and  absorbed  power  possess  sharp 
maxima  at  substantially  the  same  frequencies  as  the  maxima  for 
the  TONKS-DATTNER  resonances.  The  noise  radiation  and  the 
absorption  at  the  dipole  resonance  depend  strongly  on  collision 
frequency  and  the  plasma  tube  diameter.  Ac  low  neutral  gas 
densities  where  the  mean  free  path  of  the  electrons  is  several 
times  larger  them  the  diameter  of  the  plasma  tube,  wall  collisions 
are  more  frequent  than  volume  collisions.  Under  these  conditions 
and  if  the  diameter  of  the  plasma  tube  is  small  compared  with  the 
wavelength  (d/x  «  1),  the  shape  of  the  noise  and  scattering  spectra 
is  nearly  the  same  (see  Fig.  4).  When  the  tube  diameter  is  in¬ 
creased  the  transmission  at  the  dipole  resonance  tends  to  zeTO 
and  the  noice  radiated  decreases  rapidly.  The  shape  of  the 
spectra  for  TONKS-DATTNER  resonances,  however,  remains 
unchanged.  (See  Fig.  5. } 
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At  higher  neutral  gas  densities  die  mean  free  path 
of  electrons  becomes  smaller  than  the  diameter  of  the  plasma 
tube  and  the  volume  collisions  get  more  frequent.  The  noise 
spectrum  for  the  case  when  the  volume  collisions  are  predominant 
is  shown  in  Fig.  6.  The  noise  radiation  at  the  dipole  resonance 
increases  with  collision  frequency  and  becomes  comparable  to 
the  radiated  noise  power  at  TONKS-DATTNER  resonances. 

The  microwave  fields  in  the  waveguide  can  excite  slow 
waves  which  will  propagate  along  the  plasma  column,  and  couple 
energy  out  of  the  waveguide.  In  the  frequency  range  where  the 
dipole  resonance  occurs,  a  backward- wave  with  one  variation  of 
the  electric  field  in  azimuthal  direction  (dipole  plasma  mode) 
will  propagate.  Small  peaks  in  the  noise  and  scattering  spectra 
are  observed,  which  are  related  to  the  excitation  of  these  modes. 
(See  Figs.  5  and  6.) 


Longitudinal  polarization 

The  noise  level  for  longitudinal  polarization  is  about 

10  dB  below  the  noise  level  for  transversal  polarization.  As 

shown  in  Figs.  4,  7  and  8  the  noise  power  increases  linearly 

2 

with  the  plasma  density  (P^  -“■a)  ).  The  noise  spectrum  possesses 
maxima  at  frequencies  for  TONKS-DATTNER  resonances.  However, 
no  noise  peak  was  measured  at  the  dipole  resonance.  At  die  end  of 
the  bandpass  of  the  backward- wave  a  sharp  maximum  in  noise  and 
absorption  spectra  occurs  as  shown  in  Figs.  7  and  8. 
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VIU  MICROWAVE  EMISSION  FROM  A  CYLINDRICAL  PLASMA 


Introduction 


kT 


The  equipartition  theorem  ascribes  an  energy  to  each  degree 
of  freedom  in  a  medium  in  thermal  equilibrium.  In  a  gas  of  noninter¬ 
acting  particles  the  degrees  of  freedom  are  those  associated  with  the 
individual  particle  motion.  In  contrast  to  this  result  Bohm  and  Pines^ 
in  one  of  their  classic  papers  showed  that  in  an  electronic  plasma  some 
of  the  degrees  of  freedom  are  distributed  in  collective  modes  or  plasma 
oscillations  of  a  longitudinal  character.  As  shown  for  instance  by  the 
present  author'  '  such  longitudinal  oscillations  cannot  exist  in  bounded 
plasma  without  coupling  to  external  wave  fields.  In  other  words  one 
would  expect  thermally  excited  collective  modes  to  radiate.  The  basic 
assumption  of  the  present  work  is  that  the  radiation  reaction  has  a. 
small  effect  on  the  electronic  motion  and  that  the  emitted  power  spec¬ 
trum  can  be  calculated  assuming  thermal  equilibrium. 

Experimental  information  on  the  emission  spectrum  has  been 
published  by  Agdur  et,  al.  ^  who  measured  the  microwave  noise  radia¬ 
tion  from  a  plasma  cylinder  as  a  function  of  frequency  for  polarization 
directions  parallell  and  perpendicular  to  the  axis.  (Below  we  will 
refer  to  polarizations  as  TE  and  TM-modes  respectively  when  the 
axial  component  of  the  electric  or  magnetic  field  is  equal  to  zero. ) 

Agdur  et.  al.  measured  a  smooth  spectrum  for  scattered  TM  waves 
and  the  wellknown  resonance  structure  (see  e.  g.  Dattner^)  for  scattered 
TE  waves.  TE-noise  possessed  an  irregular  spectrum  with  several 
small  peaks  whereas  TM-noise  was  of  much  lower  intensity  and  was 
not  resolved  into  any  particular  structure.  Later  experiments  by  the 
same  group  which  have  given  a  more  detailed  information  about  these 
phenomenon  have  shown  that  .  both  polarizations  now  appear  to  have 
spectra  practically  identical  to  that  of  scattered  TE-  radiation. 


4 


V 


114. 


A  first  attempt  to  explain  these  results  was  made  in 


iJfl 


As 


to  the  scattering  spectrum  the  experimental  result s  are  now  well  ex¬ 
plained^'  The  multiple  resonances  arise  due  to  excitation  of 

longitudinal  waves  perpendicular  to  the  cylinder  axis.  The  radially 
inhomogeneous  density  distribution  of  the  electrons  strongly  modify 

the  dispersion  relation  of  these  waves  to  give  the  observed  spectrum. 

(8) 

Assuming  Rytov"s'  generalization  of  Kirchhoff  's  law  to  be  valid,  one 
should  be  able  to  use  the  same  model,  calculate  the  absorbtivity  and 
predict  the  noise  spectrum.  This  scheme  was  carried  out  in’1  ’  for 
TE-waves  emitted  perpendicularly  to  the  axis  and  lead  to  the  same 
frequency  spectrum  as  for  scattered  waves  of  the  same  kind.  As  the 
measurements  give  the  total  emitted  TE  Poynting  flux  which  contains 
contributions  from  obliquely  emitted  TS- waves,  v/e  cannot  immedia¬ 
tely  conclude  that  experimental  and  theoretical  results  agree.  It  seems 
likely,  however,  that  the  obliquely  emitted  TS- waves  will  be  coupled 
to  the  longitudinal  oscillations  and  have  approximately  the  same  fre¬ 
quency  spectrum  as  the  normally  emitted  ones.  The  purpose  of  the 
present  paper  is  to  show  what  may  be  less  evident,  that  we  will  get 
obliquely  emitted  TM-waves  with  the  same  frequency  spectrum  thereby 
explaining  the  experimental  results  for  this  polarization  direction.  The 
basic  idea  is  that  thermal  fluctuations  will  not  only  excite  collective 
oscillations  in  phase  along  the  column  but  also  oscillations  of  arbitrary 
phase  variation.  In  fact  if  the  length  of  the  cylinder  is  assumed  to  be 
infinite  the  result  of  Bohm  and  Pines  is  immediately  applicable  and  we 
should  have  equally  strong  excitation  for  all  axial  wave  numbers  up  to 
about  one  recipical  Debye  length  after  which  the  single  particle  degrees 
of  freedom  dominate.  This  argument  suggests  that  fluctuating  electric 
fields  directed  along  the  axis  will  be  generated.  The  subsequent  rather 
lengthy  calculations  will  support  this  suggestion  and  show  that  a  radia¬ 
tion  field  with  nonzero  electric  vector  along  the  axis  is  produced. 

Using  Hughyens  principle  this  result  can  easily  be  understood.  If  we 
pick  out  one  single  mode  with  axial  wave  number  of  dipole  character 
in  the  r-<p  plane  it  is  realized  with  reference  to  figure  1,  that  v/e  will 
get  a  radiation  field  of  waves  propagating  at  an  angle  3  to  the  normal 
of  the  cylinder  axis. 
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Fig.  1.,  Oblique  emission  from  dipoles  whose  phase  varies  as 

;iz.  kz  cos  Qjjj  ^  /c  ,  jj  _  2x/\  ♦  k  . 

p'  z  7  p  z 


Description  of  the  plasma 


For  our  present  purposes  it  is  adequate  to  describe  the  plasma 

by  meanB  of  hydrodynamic  equations  and  to  include  the  effect  of  the 

ions  as  a  neutralizing  background  only.  This  model  does  not  contain 

thermal  fluctuations  but  is  easy  to  handle  when  considering  geometric 

complications . 

(2) 

As  in'  '  we  will  treat  the  scattering  problem  and  then  rely  on 

Kirchhoff's  law  to  give  the  fluctuation  psectrum.  After  linearizing  and 

assuming  the  plasma  to  be  homogeneous  and  all  perturbations  to  vary 
-itot 

as  e  we  get 


1.  -ioiN  +  Nq  div  v  =  0 

2.  -iaiv  =  -  eE/m  -  (w^/Nq)  grad  N 

In  addition  we  have  Maxwell's  equations. 

3.  curl  E  =  iuJB 

4.  curl  1J  =  -  N0e"v  -  use  ^E 
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5.  div  c  QE  ■  -  eN 

6.  div  B  =  0 

Eq.  1  is  of  course  redundant  but  is  convenient  to  keep  explicitly. 

(2) 

It  was  shown  in'  7  that  Eq.  1-6  can  be  reduced  to 

7.  V2®  +(®2e/c2)B  =  0 
and 

8.  72N-H»2e/wVs0 

with 

2,2  2  /  2 

9.  =  1  -  a>“/a>  =  1  -  N^e  ye^m® 

and  appropriate  relations  between  E,  v,  B  and  N.  In  the  present 
context  we  have  found  it  more  convenient  to  introduce  the  electromag¬ 
netic  potentials  IA  and  <5.  As  usual  we  put 

10.  B  =  curl  IA 

It  turns  out  that  by  a  suitable  choice  of  gauge  the  resulting 
wave  equations  can  be  considerably  simplified.  The  key  to  this  simp¬ 
lification  is  that  as  seen  from  Eq.  4  the  current  -NQev  is  composed  of 
two  parts,  one  due  to  electron  motion  caused  by  the  average  field  and 
one  caused  by  thermal  motion  and  density  gradients.  The  first  part  of 
this  current  is  easy  to  include  in  a  frequency  dependent  permittivity, 
the  second  is  irrotational  and  can  be  removed  from  the  source  term  by 
a  suitable  gauge  transformation.  Thus  if  we  put 

11.  4>  =  -fv2e/s2ECg)  N 

it  is  easy  to  show  that  the  resulting  wave  equations  are 

12.  ^lA-Ko^e/c2)  IA  =  0 

13.  +(®2c/w2  )  9  =  0 

14.  div  IA  =  0 

15.  E  =  -kuIA  -  7* 

The  transversality  condition  14  shows  that  eq.  12  represents 
transverse  waves  which  in  the  absence  of  boundaries  are  uncoupled  to 
the  longitudinal  waves  represented  by  Eq.  13.  Boundary  conditions 
are  as  usual  continuity  of  the  tangential  components  of  the  E  and  B 
fields.  In  addition  because  of  the  presence  of  a  longitudinal  component 
we  need  a  boundary  condition  on  9.  This  is  obtained  from  die  require¬ 
ment  that  the  normal  component  of  the  particle  current  must  vanish  at 
the  boundary.  Expressing  this  current  in  terms  of  the  potentials  we  get 

16.  7*  -t  q(®21A  +  ioV*) 
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Oblique  scattering  of  TM- waves 


For  obli<pieincidencs  of  a  TM  wave  we  will  get  excitation  of 
both  TM  and  TE  components.  Inordex  to  derive  suitable  exprasnionn 
for  the  TE  waves  we  introduce  the  superpotential  jTdefined  by 


17. 


t.TE 

IA.  s 


curl  u 


In  addition  to  Eq:s  12  -  15,  describing  the  fields  in  die  plasma, 
we  have  die  usual  free  space  wave  equation  outside  the  cylinder.  De¬ 
composing  the  solutions  into  TM  and  TE  components  and  choosing 
boundary  conditions  at  infinity  corresponding  to  an  incoming  plane 
wave  and  an  outgoing  scattered  wave,  we  get  in  vacuum: 


18.  y  A 


ik  i(nq>  +  k  z) 

Ar  *  -“T-  5  IJn(a  r)  +  an  Hn(a r)  1  e 


TM 


,  TM 


k  ,  .  i(n*P  +  k  *) 

— 5—  E  n  [ J  (a  r)  +  a  H  (a  r)  ]  e 
2n  n'  '  nnv 
a  r 


i(n<p  +  k  z) 

A*"*  -  Z  [j  (ar)  +  a  H  (ar)]  e 
z  n  nv  '  n  nx  ' 


TM 
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TE  ik  A 

ix™  =  — 2L_  z  b  H’(ar)e  * 

n  n  nv  9 


TM  kza  „  .  „ ,,  , 

~ - r  &  n  b  H ar  e 

2  n  n  n 
a  r 


i(nq>  +  kzz) 
i(n<J>  +  k^z) 
i(n?  +  k  z) 


u  ™  =  a  E  b  H  (a  r)  e 
r  z  n  n  n' 


a  denotes  the  cylinder  radius  and  we  have  also  introduced  the  abbre¬ 
viations 


20.  a  =  —  cos  9 ,  k  =  ^  sin  9 
c  z  c 

and  dropped  the  superscript^  on  the  Hankel  functions,  a^  and  bQ  are 
arbitrary  constants  to  be  determined  by  die  boundary  conditions.  Prime 
denotes  derivative  with  respect  to  argument.  To  get  the  transverse 
fields  in  the  plasma  we  use  the  same  expression  with  die  modification 
that  new  constants  (d  and  f  )  are  introduced,  die  Hankel  junctions  are 
replaced  by  Bessel  functions  regular  at  the  origin,  the  incoming  wave 
term  is  omitted  and  a  in  die  argument  is  replaced  by  o  *  defined  by 

,2  /  .  2  _ . 

a  =  — ^  (*  -  sin  9) 

c 
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We  also  get  for  the  scalar  potential 

i(nq>  +  k  z) 

22.  •  «  2  gn  Jn(  pr)  e  z 

with  2  2 

23.  p2  -  ^2  (*  -  sin20) 

wfc  c 

Applying  boundary  conditions  at  r  *  a  we  get  die  following  set 

of  equations  for  the  determination  of  the  constants  a  ,  b  ,  d  ,  f  and 

&  n  n.  n 


^  [  J  (a a)  +  a  H  (a a)]-^b  H'(a a)  -  J  (•  ’ a)d- 

*  **  n  n  n  ■£,  n  n  . n  n 

a  a  ac  a  a 


2 

a)  a  c 


r(a'a)fn-  n  Jn(Pa)  gn 


J  (a  a)  +  a  H  (a  a)  =  J  (a  *  a)  d  -  ak  J  (Pa)  g 
n'  7  n  n'  7  nv  7  n  z  n'^  7  6n 


k  n  , 

_*  u  .i  v.  *  r 


'  TT  Hn{aa>bn  '  «  [j;(aa)  +  an  Hk(aa)1  *  ’  ^  dn  ‘ 


c  n  J  (a  ’  a)  f 
n  7  n 


z  J’(u  ’  a)  d  +  -  j  2n  J  (a  ’  a)  f  +  <D2pag  J’(pa)  =  0 
n  n  i  a  n  n  n  n 


2  2 
(II  to  c 


b  H  (o  a)  =  t  f  J  (a  *  a) 
n  nx  7  n  n'  7 


Solving  for  b^,  the  amplitude  of  the  scattered  TE-waves,  we  get 


25.  b  *. 
n 


Jn(aa)Qn 


(kza)  Hn(aa)  ^ Pn  +  TT  *  ~TT^  *  ?!  1 

n  kz 

where 

H’(aa)  J'(oa) 

■n.  n  _  r  nx  7  n'  7  i  a  a 

26,  Qn  -  lH"(oaJ - T  (aafJ  ~ 

n'  n'  7 

,  x2  H’  (oa)2  2  ,  2  J' (a 7  a)H7  (o a) 

„  n  .  .  «»  n'  7  a  r,  aa  a  n'  7  n'  7  t 

27‘  Pn  ~  1 - T~  —  -^-—72  11 - Z - H  (oa)J  (o' a)  J 

n  H  (aa)  a’  n  n'  7  n'  7 


28.  S 


gar)!..)  t  r (. ■ a)  a,t  h;(..)  -V  , 

n  nVPa*  SJa '  *}  n  T*"2  n  Hn<oa*  c  (o'a)2 
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j;<a'a) 


ta 

n  TTfaXI 


30. 


j;(a'a)  paj;(pa) 
T^aJ-  nr(W 


To  find  possible  resonance  conditions  we  now  study  die  expres¬ 
sion  25  in  the  limit  <ua/c  «  1  i.  e.  when  the  vacuum  wavelength  be¬ 
comes  much  larger  than  the  plasma  radius.  From  the  Taylor  expansions 
of  the  Bessel  functions  involved  it  is  easily  seen  that  unless  the  real  part 

of  the  bracket  in  the  denomination  of  25  vanishes,  b  is  of  the  order 

2n  —  1  ® 

(coa/c)  .  (We  disregard  the  factor  (k  a)~  which  originates  from  the 

z 

fact  that  at>n  is  the  amplitude  of  ]I  whereas  the  amplitude  of  IA  in  the 
incoming  wave  is  put  equal  to  unity)  If,  however. 


nJn(Pa)  2m 

s-r 


the  real  part  of  the  bracket  vanishes  and  b  becomes  of  the  order  of 

oa  n 

unity.  We  observe  that  in  the  limit  —  «  1  no  dependence  on  die  angle 
of  incidence  appears  in  equation  31.  The  resonance,  condition  31  coin¬ 
cides  with  the  resonance  condition  obtained  in^  for  perpendicularly 
incident  TE-waves.  This  should  not  be  surprising  as  the  resonances 
are  caused  by  longitudinal  waves  perpendicular  to  the  cylinder  axis. 
These  waves  are  coupled  to  transverse  fields  only  by  a  radial  density 
discontinuity  at  the  boundary.  Keeping  this  physical  picture  in  mind, 
one  would  not  expect  a  phase  difference  along  ?he  cylinder  axis  to  alter 
the  resonance  conditions.  This  lack  of  dependence  on  the  angle  of 
incidence  is  important  to  remember  when  discussing  experimental 
results.  As  we  will  outline  below  the  noise  frequency  spectrum  will 
have  peaks  when  31  is  satisfied.  If  there  were  a  strong  angular  de¬ 
pendence,  this  would  show  up  as  a  broadening  of  the  peaks  in  the  mea¬ 
sured  spectrum  which  is  the  cos  6 -weighted  sum  of  all  waves  emitted 
with  a  component  of  die  electric  vector  along  the  cylinder  axis.  In  the 
experiments  of  Kerzar  there  is  no  evidence  of  any  additional  broading 
in  the  noise  spectra  as  compared  to  the  scattering  spectrum  for  nor¬ 
mally  incident  TE-waves. 
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It  should  also  be  noted  that  for  normally  incident  TM-waves 
2  2 

the  factor  a  /k  in  die  denominator  of  Eq.  25  will  make  b  equal  to 
z  n 

zero  even  at  frequencies  when  31  is  satisfied.  For  completeness  we 
give  the  scattered  TM  amplitude  for  normal  incidence 


Jn(a  a) 


hr  3  tail  T  (a T a]  1 

_ n _ n _ 

,  J'(a’a)  _  H'(oa) 

ra  a  n'  '  nan'  i 

^  n  J"  (a  ’  a)  ’  n  It  (a  a)  ' 
n  n 


which  is  seen  to  possess  no  resonances  and  tend  rapidly  to  zero  with 

toa/c. 

Noise  spectrum 

Assuming  Rylov's  result  to  apply  ^  we  should  now  be  able  to 
calculate  the  spontaneously  emitted  energy  per  unit  surface  as  frequency 


d  kT  A 
P  ~~ — 2 A“ 
®  2rX 4 


a  -  /  Re 
cu  -  /  Re 


£  dv 


where  index  i  stands  for  incident  radiation. 


Now  it 


shown  in'  ’  that  for  the  present  model  we  have 


r  *  c  u\c 

j  Re^*  E  dv  =  j 


*)  This  is  an  assumption  to  be  verified  by  die  results  as  Kirch- 
hoff's  law  is  derived  either  in  die  approximation  of  geometrical  optics 
or  by  Rytov  in  media  where  the  relation  between  £  and  1)  is  local  and 
linear.  Our  model  is  equivalent  to  a  nonlocal  linear  relation. 


V 
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The  collision  frequency,  v  ,  is  defined  as  the  recipical  mo- 

c 

mentum  relaxation  time.  A  term  accounting  for  this  dissipation  must 
now  be  added  to  equation  2.  Inserting  16  in  35  and  using  same  simple 
vector  identities  we  get 


37. 


A©=* 

2 

to  e , 


2  4 
c  -  w 

JL_E_ 


r 

j 


ia. 


dv  + 


®4,«o 


w 


r 


|  *  | 2  dv  + 


.  2  2 

r  *  i*  »©  r  +  * 

j  «  V4>  *  dl>  +  -P  j  (IA  «  -  IA*  )  d? 


(2) 

As  was  done  in'  we  now  assume  that  IA  and  *  can  be  obtained 
neglecting  v  It  is  then  a  straight  forward,  but  tedious,  procedure 
to  insert  IA  and  $  from  17,  18  and  22,  carry  out  all  integrations  which 
turn  out  to  be  elementary  and  insert  constants  of  integration  from  24. 

We  omit  these  calculations  here  and  quote  the  main  result  only.  The 
absorbtivity  possesses  sharp  maxima  at  frequencies  given  by  31, 
if  ©a/c  «  1.  Consequently  we  have  shown  that  in  thermal  equilibrium 
the  noise  radiation  with  the  E-vector  parallell  to  the  cylinder  axes 
should  possess  the  same  frequency  spectrum  as  for  instance  scattered 
TE-waves.  There  also  seems  to  be  good  reason  to  believe  that  the 
noise  spectrum,  as  has  been  shown  for  die  scattering  spectrum,  should 
be  strongly  modified  by  the  non-homogeneous  radial  density distr.  and  ten 
agree  well  with  experimental  results.  The  measured  noise  spectra 
therefore  appears  to  give  good  evidence  for  the  result  of  Bohm  and 
Pines  that  part  of  the  kinetic  energy  in  a  plasma  is  distributed  over 
the  collective  degrees  of  freedom. 
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